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PREFACE 


This  final  technical  report  covers  a study  program  conducted  to  define 
important  flight  control  system  design  and  handling  quality  criteria  for 
moving  loads  slung  beneath  tandem-rotor  heavy  lift  helicopters.  The 
study  was  conducted  by  members  of  the  Flight  Simulation  Test  Group, 
with  Mr.  J.B.  Sinacori  serving  as  Principal  Investigator.  Principal 
contributors  to  the  simulation  portion  of  the  study  were  Messrs. 

A.Y.  Murakoshi,  R.  Silvestri,  W.G.  Spring,  and  W.L.  Ross.  The  non- 
linear sling  loads  analysis  was  performed  by  G.L.  Danek. 

The  basic  study  involved  flight  simulation  of  the  Model  347  advanced 
tandem-rotor  helicopter  and  was  sponsored  by  the  U.S.  Army  Air  Mobility 
Research  and  Development  Laboratories  (USAAMRDL),  under  Contract  DAAJ02- 
72-C-0112.  Mr.  Robert  P.  Smith  served  as  USAAMRDL  monitor.  Work  was 
started  in  July  1972  and  the  study  was  completed  in  October  1973.  During 
that  15-month  period,  a detailed  mechanization  of  a tandem-rotor  heavy 
lift  helicopter  was  developed  for  use  with  Northrop' s large-amplitude 
moving-base  simulator  and  its  integrated  computer  facilities.  These 
simulation  facilities  and  the  detailed  helicopter  mechanization  were  used 
during  the  15-month  contract  period  on  a directly  related  USAAMRDL  study 
(Contract  DAAJ02-72-C-0047)  involving  the  stabilization  of  externally 
slung  helicopter  loads  and  for  related  HLH  control  system  studies 
conducted  by  Boeing -Vertol  and  a closely  related  study  conducted  by 
Northrop  under  direct  contract  with  Boeing. 

As  work  on  these  closely  related  study  projects  progressed,  the  common- 
ality of  many  study  objectives  became  apparent.  The  impact  of  this  common- 
ality resulted  in  modification  of  certain  aspects  of  the  Contract  DAAJ02- 
72-C-0112  study  to  avoid  duplication  of  effort  and  to  take  advantage  of 
the  detailed  simulator  mechanization  set  up  for  the  other  directly  related 
control  system  simulation  studies.  Thus,  a substantial  body  of  data 
developed  dui ing  the  conduct  of  the  Contract  DAAJ02- /2-C-0047  study  and 
certain  data  developed  during  the  Boeing-Vertol  simulator  studies 
have  been  integrated  with  Contract  DAAJ02-72-C-0112  data  and  are  contained 
in  this  report. 

The  major  subsections  that  comprise  the  body  of  this  report  discuss 
derivation  of  the  equations  of  motion,  analysis  of  sling  load  carrying 
problems,  the  dynamics  of  an  8-by-8-by- 20-foot  container  slung  beneath 
a tandem-rotor  helicopter,  and  a Model  347  helicopter  flight  simulation 
with  external  loads.  Methodologies  used  in  developing  this  data  include 
theoretical  analyses,  acquisition  and  evaluation  of  static  and  dynamic 
wind  tunnel  and  flight  test  data,  and  piloted  ground-based  simulation. 

The  discussion  is  followed  by  a brief  statement  of  conclusions  reached,  and 
the  report  concludes  with  an  appendix  containing  functional  block  diagrams 
and  other  data  relative  to  the  computerized  flight  simulation  of  the 
Model  347  helicopter. 
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AIRFRAME  DYNAMIC  MODEL  AND  EQUATIONS  OF  MOTION 


The  mathematical  model  used  for  simulation  of  the  full  flight  envelope  of 
the  Model  347  tandem-rotor  heavy  lift  helicop'.er  is  described  in  detail 
in  Reference  1 and  is  summarized  in  this  section.  Brief  descriptions  of 
each  major  calculation  are  presented  in  the  general  order  of  their  per- 
formance, starting  with  equations  of  motion  and  concluding  with  the 
forces  and  moments  that  drive  the  ..totion  equations. 

The  helicopter  is  modeled  as  a rigid  body,  with  symmetry  assumed  about 

the  X-Z  plane  so  that  I and  I are  zero.  Linear  and  angular 
r xy  yz 

accelerations  of  the  airframe  are  expressed  in  the  helicopter  body  axis 
system.  The  forces  X,  Y,  and  Z and  moments  L,  M,  and  N contain  only 
aerodynamic  terms.  Euler  angle  rotations  are  performed  to  orient  the 
airframe  in  inertial  space.  The  basic  equations  of  motion  providing  linear 
and  angular  accelerations  about  the  X,  Y,  and  Z axes  are: 

Linear  Acceleration 

u = ^ - g sin  0 + Rv  - Qw 
m 

to  = ^ + g sin  A cos  0 + Pvv  - Ru 
m ' 

to  = “ + g cos  A cos  0 + Qu  - Vv 
m ' 


Anr^ilar-Aceelera^ijan. 


1.  Cogati,  C.,  Gajkowski,  B.J.,  and  Gumett,  T.S.,Jr. , FULL  FLIGHT 

ENVELOPE  MATHEMATICAL  MODEL  FOR  347/HLH  CONTROL  SYSTEM  ANALYSIS,  The 
Boeing  Company  Report  No.  D301-10148-1,  December  1972. 
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The  linear  and  angular  accelerations  of  the  helicopter  were  mechanized 
for  the  liybrid  simulation  on  an  analog  computer.  By  integrating  these 
accelerations  with  time,  linear  and  angular  velocities  in  the  body  axis 
are  produced.  Fuselage  angle  of  attack  and  sideslip  information  are 
derived  from  body  axis  linear  velocity  resolutions,  corrected  for  the 
influence  of  rotor  downwash  on  the  fuselage  and  for  steady  wind  or  gust 
velocity. 

The  total  linear  aerodynamic  velocities  associated  with  fuselage  aFy<j  and 
3 consist  of  body  axis  ground  speed,  steady  wind  velocity  resolved 
into  the  aircraft  body  axis,  and  gust  velocity  resolved  into  the  aircraft 
body  axis.  Steady  wind  velocities  introduced  into  the  linear  velocity 
equations  are  based  on  a steady  wind  (V  ) blowing  from  a fixed  azimuth 
(Y  ) in  the  inertial  X-Y  plane  on  to  the  helicopter.  This  wind  from  the 
fixed  inertial  axis  is  resolved  into  the  body  axis  and  rotor  axes  of 
the  helicopter. 

Before  the  fuselage  a„gle  of  attack  can  be  defined,  using  aerodynamic 
velocities,  the  effect  of  rotor  downwash  on  the  fuselage  must  also  be 
considered.  The  forward  and  aft  rotors  both  produce  a downwash 
component  on  the  fuselage  that  is  assumed  to  be  equal  to  an  ’’average" 
induced  velocity  of  the  rotor  system.  Fuselage  sideslip  angles  do  not 
require  correction  for  rotor  downwash  interference  and  may  therefore  be 
defined  as  functions  of  total  longitudinal  and  lateral  aerodynamic 
velocity. 

Rotor  forces  and  moments  can  be  calculated  as  soon  as  the  rotor  control 
inputs,  angle  of  attack  and  sideslip  values,  and  both  inflow  and  advance 
ratios  are  known.  Rotor  angle  of  attack  and  sideslip  are  computed  from 
linear  velocity  components  in  the  same  way  as  described  for  the  fuselage. 

The  first  step  in  computing  rotor  sideslip  or  angle  of  attack  is  to  sum 
the  total  velocity  components  with  the  velocity  of  each  rotor  relative 
to  the  aircraft  center  of  gravity  created  by  the  aircraft's  angular 
pitching  motion  (tangential  velocity  component  resulting  from  airframe 
angular  velocities). 

Next,  the  linear  velocities  at  each  rotor  hub  are  resolved  through  the 
shaft  incidence  angle  into  the  shaft  norma]  plane  (SNP).  According 
to  classical  rotor  theory,  all  rotor  force,  moment,  and  flapping 
calculations  are  based  on  the  assumption  thac  the  rotor  faces  directly  into 
the  wind,  with  no  sideward  velocity  component  existing.  To  meet  this 
requirement  for  zero  sideward  velocity,  the  longitudinal  and  lateral 
velocity  components  at  the  rotor  hubs  are  resolved  through  the  rotor  side- 
slip angle  (3')  into  the  SNP  wind  axis.  A single  velocity  vector  in  the 
rotor  SNP  is  formed  and  is  defined  as  being  parallel  to  the  local  wind 
direction  in  the  rotor. 
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Rotor  advance  and  inflow  ratios  and  \)  are  defined  as  the  ratio  of 
the  rotor  SNP  wind  axis  longitudinal  and  vertical  velocity  components, 
respectively,  to  the  rotor  tip  speed.  These  values  of  advance  ratio  are 
used  for  all  rotor  calculations  carried  out  with  the  classical  equation 
approach.  Rotor  maps,  on  the  other  hand,  require  the  use  of  a modified 
advance  ratio  because  the  tabular  map  data  is  nondimens ional. 

Inflow  ratio  is  comprised  of  terms  resulting  from  free-stream  velocity 
(normal  to  the  rotor)  and  induced  velocity.  Two  terms  are  added  to  the 
free-stream  velocity  component  of  the  inflow  ratio.  One  accounts  for 
the  uniform  induced  velocity  of  the  rotor  under  consideration, and  the 
second  reflects  the  interference  downwash  of  the  other  rotor  imposed  on 
the  rotor  under  consideration. 

Each  rotor  of  a tandem-rotor  helicopter  receives  only  three  types  of 
input  under  pilot  control! 

. Collective  Pitch  0 , produced  by  longitudinal  and  collective  stick 

or  Stability  Augmentation  System  (SAS)/Automatic  Flight  Control 
System  (AFCS)  T'.puwS. 

. Longitudina’.  Cyclic  (B  „),  normally  comprised  of  "q"  sensed  cyclic 
trim  and/or  AFCS  cyclic  "on  the  stick"  inputs. 

. Lateral  Cyclic  (A^q),  comprising  a combination  of  lateral  stick 
and  pedal  or  AFCS/SAS  inputs. 

Rotor  input  collective  and  cyclic  controls  (9q,  B^,  A^)  are  compiled 

and  both  longitudinal  and  lateral  cyclics  are  resolved  through  the  rotor 
sideslip  angle  to  align  controls  orthogonally  to  the  local  wind.  Forward 
rotor  cyclic  inputs  are  corrected  for  delta  three  (6^)  hinging  effects. 

These  corrections  account  for  the  phasing  of  the  control  system  to 
produce  orthogonal  flap  response  and  for  the  change  in  blade  pitch 
occurring  because  of  the  6^  pitch-flap  coupling  effects. 

Coning  (aQ),  longitudinal  (a^),  and  lateral  (h^)  first  harmonic  flapping 

coefficients  are  computed  and  first-order  dynamic  lags  are  applied  to 
account  for  rotor  system  flapping  response.  Representation  of  rotor 
flapping  in  the  dynamic  model  is  accomplished  by  using  quasi  steady-state 
classical  calculations  based  on  the  following  assumptions: 


First  harmonic  theory 
Uniform  inflow 
No  reverse-flow  effects 
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. Identical  forms  for  front  and  rear  equations 
. Uniformly  twisted  untapered  blades 

. No  6,  effects  (accounted  for  by  correcting  input  controls, 
as  discussed  earlier) 

. Zero  tip-loss  factor 
. Negligible  hinge  offset 
. Rigid  blades 

. No  stall  and  compressibility  effects 

. Constant  lift-slope  and  drag  coefficients  varying  only  with  C^. 

Rotor  forces  and  moments  are  computed,  using  either  the  classical  equation 
method  or  the  rotor  map  approach.  In  the  former  case,  forces  are  resolved 
perpendicular  and  parallel  to  the  SNP  and  aligned  with  the  SNP  wind  axis. 
In  the  latter  case,  rotor  thrust,  power,  and  propulsive  forces  are 
derived  from  isolated  tabular  rotor  or  wind  tunnel  data  and  expressed  in 
the  rotor  wind  axis. 

Fuselage  aerodynamic  forces  and  moments  ere  calculated  from  tabular  data 
derived  from  wind  tunnel  ttsts.  Rotor  aerodynamic  forces  and  moments  are 
resolved  through  rotor  sideslip  and  shaft  incidence  angles  back  into  the 
aircraft  body  axis.  These  resolved  forces  and  moments  are  summed  with 
body  axis  (fuselage)  forces  and  added  to  each  moment  imposed  on  the  air- 
craft by  individual  rotor  T,  H,  and  Y forces.  These  force  and  moment 
sums  are  inserted  in  the  equations  of  motion  to  compute  the  airframe 
acceleration  set. 

Engine  governor  response  to  aerodynamic  rotor  power  requirements  is 
computed  along  with  individual  rotor  speed  degrees  of  freedom.  Rotor 
shaft  windup  spring  rates  are  accounted  for  in  the  governor  model. 

The  external  sling  load  is  coupled  with  the  rigid-body  helicopter 
equations  of  motion  described  on  page  6.  The  aerodynamic  force  and 
moment  data  used  with  the  sling  load  equations  is  given  in  Ta> le  1. 
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(Sharp-Edged  Box  - Body  Axis) 

YAWING  MOMENT,  N^  SIDE  FORCE,  YLT 


0 |54.94  62.99  67.72 

57.27  63.69  66.70  42.7 

58.76  64.39  69.68  64.8  

60.73  66.62  70.47  88.3  -27.0  -89.4  I -21.6 
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-59.4 

-67.0 

-82.1 

-99.0 

-115.6 

-144.1 

-180.9 

-225.7 


-4.0 

-19.8 

-35.5 

-49.9 

-61.5 

-68.4 

-75.8 

-85.9 

-93.4 

-103.5 

-113.0 

-122.6 

-155.3 

-184.8 

-225.7 


NOTE:  Units  for  the  aerodynamic  forces  and  moments  above  are: 

Dlt(FT2)  = D/q,  Nlt(FT3)  = N/q,  and  Ylt(FT2)  = Y/q. 


The  method  of  integrating  the  sling  load  equations  with  those  of  the 
helicopter  is  described  in  following  paragraphs. 


' .V  V~-  ■■■  . ~ •'  : 


The  rigid-body  equations  for  two  vehicles  (helicopter  and  load)  were 
linked  by  additional  equations  describing  cable  tension  forces.  These 
forces  and  their  corresponding  moments  were  added  to  the  forces  and  moments 
due  to  gravitational  and  aerodynamic  effects.  The  key  is  to  describe 
the  movements  of  the  lower  attach  point  relative  to  the  upper  one  for 
each  cable,  and  compute  the  tension  force  as  a function  of  the  cable 
elongation  heyono  the  zero  tension  value.  The  cable  forces  are  written 
in  terms  of  the  stretched  length  and  length  rate  and  are  used  together 
with  matrix  representations  for  the  displacements  of  the  lower  attach 
points  (load)  relative  to  helicopter-fixed  axes.  The  coupled  equations 
describing  the  movements  of  two  vehicles  in  twelve  degrees  of  freedom 
when  there  is  a single  cable  attached  to  both  vehicles  follow: 


The  nomenclature  describing  the  axes  systems  used  and  attach  points  for 
one  cable  is  shown  above.  In  this  diagram,  l is  the  straight  line 
distance  between  the  attach  points.  When  this  distance  l is  greater 
than  L , the  zero  tension  cable  length,  a tension  T exists  which  will 
produce  forces  and  moments  on  helicopter  and  load, thereby  coupling  them 
together  dynamically.  It  now  remains  to  compute  the  angles  and 
and  postulate  the  cable  tension  T as  a function  of  jf  and  its  derivatives. 
The  format  for  these  computations  shown  for  nth  number  of  cable  systems 
is  given  on  the  following  pages. 
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1)  Difference  between  lower  attach  point  velocity  and  upper  attach  point 
velocity,  relative  to  helicopter  axes  for  the  nch  cable. 


UBA 

UB  + VB  * Vb 

“a  + Va  • Va  " 

n 

VBA 

n 

■ m 

n n 

VB  + rBXB  ' VB 
n n 

- 

n n 

VA  + Va  ■ Va 

n n 

WBA 

WB  + PByB  “ qBXB 

“A  + pAyA  • qAXA 

n 

n n 

n n 

2)  Difference  between  attach  points,  computation  of  j^,  X^n»  X2n  an^ 
their  rates. 


I 


( = JUBA  dt  * YBA  = ^VBA  dt  , Z 

BA  o n n o n ’ BA 

f {BA  \ /"YBA 


* Kdt 


4„  * «BA  + YBA  + 2BA 
n n n 

V /-i  ; \ 

■ C +\ " + 2nj 


tan  Au  tan  X2n 


-Y  * 

. BA  j l 

X2n  = V°sX2n  ‘ ti  t3n  X2n 


^n  XBA  sinXln  CosX2n  " YBA  sin?‘2n 
n n 


+ ZBA  “s*ln  cos*2n 
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3)  Resolution  of  the  nth  cable  tension  into  forces  and  moments  in 
helicopter  and  load  axes. 


Y*r  * [»]  [a]T  [\,]T  0 , "ox  - yb  x ybt 

nLJLJLJ  » n n n 


B BT 

n n 

_J  L.  - 


yA  y YaT 
n X n 


In  the  above  expressions,  the  matrices  [a]  and  b]  represent  the  commonly 
used  orientation  matrices  of  direction  cosines  jor  the  helicopter  ( [aJ  ) 
and  the  load  ( ) respectively.  The  matrix  \n  represents  the 

transpose  of  the  orthogonal  matrix  of  direction  cosines  of  the  taut  cable 
after  being  rotated  through  the  angles  and  \2n» 


■>  in\  2n 


-sim.  cosX,  sin?u  cosX.  cosX,, 


4)  Finally,  the  «,able  tension  equations  must  be  written.  In  the  case 
chosen,  the  cable  is  assumed  to  tense  when  its  length  is  greater 
than  its  zero-tension  value.  The  tension  functions  are  assumed  to  be 
linear;  i.e.,  the  force  is  proportional  to  the  cable  stretch  and 
stretch  rate. 


(CABLE  SLACK)  (CABLE  TAUT) 


Other  tension  functions  may  be  used  to  represent  sling  material  such  as 
nylon, which  is  known  to  have  nonlinear  tensile  properties. 
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ANALYSIS  OF  SLING  LOAD  CARRYING  PROBLEMS 


The  validity  of  the  linear  analysis  of  sling  loads  was  tested  by  comparing 
its  results  with  flight  data.  A number  of  candidate  sling-load  flight 
tests  were  examined,  and  two  well -documented  cases  were  selected  for 
comparison.  The  first  case  involved  a test  of  an  empty  MILVAN  container 
on  a two-cable  suspension  by  a CH-47C  helicopter  as  described  in 
Reference  2.  The  container  was  carried  small-end-into-the-wind,  using 
a spreader  beam  weighing  3505  pounds.  Reports  of  the  pilot  and  observer 
about  container  motions,  and  the  pilot's  judgments  about  maximum  allowable 
carrying  speed,were  studied.  The  results  are  summarized  in  Table  2. 

The  second  case  involved  a lift  made  in  Colorado  by  a CH-54  Tarhe. 

The  load  was  a 13,500-pound  mobile  home  slung  on  a single  cable  with 
cross-slings  and  spreader  bars.  Observer  reports  from  a chase  helicopter 
and  film  records  were  examined  in  detail.  The  findings  are  summarized 
in  Table  3. 

In  general,  the  comparisons  are  very  poor,  particularly  when  the  container 
is  carried  small-end-forward  on  a two-point  suspension.  The  linearized 
analysis  predicted  no  instabilities,  with  an  increasing  yaw  frequency  as 
airspeed  increases,  while  the  observations  showed  sustained  yaw  oscillations 
t about  the  zero  airspeed  (bifilar)  frequency.  In  the  broadside  case, 

.iO  pilot  reports  of  oscillatory  frequency  were  available.  Indeed,  the 
analysis  predicted  a damped  oscillatory  motion  which  was  not  apparent 
in  the  film  records.  The  comparison,  therefore,  is  quantitatively  incon- 
clusive but  qualitatively  interesting  because  t.»e  load  did  stabilize 
to  a broadside  position  and  the  airspeeds  flown  were  less  than  the 
predicted  maximum  stable  airspeed. 

This  generally  poor  correlation  coupled  with  the  lack  of  good  data  for 
dynamic  comparison  prompted  the  setup  of  an  analog  computer  to  conduct 
a nonlinear  analysis.  A parallel  theoretical  effort  was  also  initiated. 

The  simple  linearized  equations  of  motion  previously  described  were 
programmed,  the  yaw  axis  was  made  to  allow  continuous  rotation,  and  yawing 


2.  Watkins,  T.C.,  Sinacori;  J.B.,  and  Kesler,  D.F.,  STABILIZATION  OF 
EXTERNALLY  SLUNG  HELICOPTER  LOADS,  Northrop  Corporation  Aircraft 
Division;  USAAMRUL  Technical  Report  73-42,  Eustis  Directorate, 

U,  S,  Army  Air  Mobility  Research  and  Development  Laboratory,  Fort 
Eustis,  Va.,  August  1974. 
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TABLE  2.  FLIGHT  DATA  ANALYSIS  COMPARISON 
(8-By-8-By-20-Foot  Empty  Container,  Two-Point  Suspension) 


FLIGHT  DATA 


ANALYSIS. 


EO-FT 

CABLES 


8-FT 

SLINGS 


CH-47C 


SPREADER 

8-X-8-X-20-FT 
CONTAINER 


Load  Weight 

Container 5200  lb 

Spreader  Beam 3505  lb 

Cable  Length 50  ft 

Sling  Length 8 ft 

Cable  Separation 24  ft 

Rigged  Pitch  Angle 0° 

OBSERVATIONS 

The  level  container  with  50-ft 
cables  exhibited  a yaw  oscillation 
of  0.25  Hz  (1.57  rad/sec)  with  an 
amplitude  of  ±35°  at  an  airspeed 
of  90  kn.  This  load  oscillation 
coupled  with  a lateral  swinging 
induced  an  aircraft  roll  oscilla- 
tion of  ±5°.  The  load  motion  fed 
forces  into  the  aircraft  which  were 
heavy  and  uncomfortable.  The 
maximum  recommended  airspeed  for 
operational  conditions  was  con- 
sidered to  be  70  kn. 


Load  Weight 87Q5  lb  „ 

Estimated  Yaw  Inertia  .9011  slug-ft^ 

l 62  ft 

50  ft 

XC 12  ft 

Angle  of  Attack 0° 

C = + 1.74/rad  (-10°  <3  < 10°) 

ne 

C = -2.06/rad  (-10°  < 3 < 10°) 


C = -1.625 
n 

r 


(-10°  < 3 < io°) 


ROOT  PLOT 

ju; 


FOR  1YAW 

INCREASING 
AIRSPEED 

'sway 


D 0.721  V 


No  instability  occurs  up  to  airspeeds 
of  100  kn. 

*Estimated  from  the  assumption  of  a 
longitudinally  uniform  container 

.2 

iz  L , x =io  ft 

C 3 BL 


SYNOPSIS 


Results  of  the  analysis  show  no  instability  with  damped  yaw  oscillation  of 
frequency  that  increases  with  airspeed.  Sway  damping  increases,  but  fre- 
quency is  nearly  constant.  Flight  observation  shows  a sustained  yaw 
oscillation  at  the  zero  airspeed  frequency  *of  the  analysis.  The  yaw 
amplitude  is  obviously  exceeding  the  sideslip  constraint  for  the  aero- 
dynamic data  of  the  analysis.  Generally,  there  is  very  poor  correlation 
between  analytical  results  and  flight  observations. 


TABLE  3.  FLIGHT  DATA  ANALYSIS  COMPARISON 
(Mobile  Home,  Single-Point  Suspension) 


FLIGHT  DATA 

£T 

»CH-54 

__  CABLE 

SPREADER 
SLINGS 

MOTOR  HOME 

Load  Weight 13,500  lb 

Load  Dimensions 

Length 23  ft 

Width 7 ft 

Height 8 ft 

Cable  Length 

Hook  to  load  ctr  of  mass...  40  ft 
Rigged  Pitch  Angle -5° 

OBSERVATIONS 

The  load  tended  to  rotate  in  yaw, 
but  stabilized  in  the  broadside 
position.  Maximum  airspeed  reached 
was  about  60  kn.  No  instability 
was  encountered,  even  in  turbulence. 


ANALYSIS 

Load  Weight 13.500  lb 

Estimated  Yaw  Inertia  18,482  slu£-ft^ 

40  ft 

X 0 

c 

Angle  of  Attack -5° 

C = 0.87/rad  (±80°  < g.<  ±100°) 

no  L 

CYP  = -0.07/re  1 (±80°  < pL<  ±100°) 

C P = -1.625  (±80°  < pL<±100°) 

**  y2 

'VS 

*h  ■ ~ r1  ’ \ - u-5  ft 

Li 


ROOT  PLOT 


JU) 


-V  - M , f- 

KNOTS  I = 0, 


SWAYW  0.897 


I YAW 


SWAY  MODE  BECOMES 
UNSTABLE  AT 


68  KNOTS 


&YKQESIS. 

The  significance  of  this  comparison  is  that  the  load  stabilized  in  the 
broadside  position,  which  is  verified  by  the  analysis  in  which  two  broad- 
side yawing  moment  nulls  are  predicted,  with  dynamic  stability  about 
those  points  below  an  airspeed  of  68  knots.  A specific  yaw  frequency 
check  is  not  possible  because  it  was  not  observed  in  flight.  This 
case  indicates  positively  that  the  static  aerodynamics  used  with 
a simple  analysis  may  in  fact  be  predicting  reality  correctly. 
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moment  data  for  any  value  of  yaw  and  sideslip  angle  was  incorporated. 

This  extension  of  the  analysis  did  not  provide  results  that  compared 
favorably  with  flight  data.  No  sustained  oscillation  could  be  produced, 
although  many  of  the  properties  of  the  broadside  case  could  be  reproduced. 
Therefore,  a redirection  of  effort  became  necessary  and  is  described  in 
the  following  section. 
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DYNAMIC  WIND  TUNNEL  TESTS 
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TEST  PROCEDURE 

Previous  attempts  to  predict  the  action  of  sling  loads  by  using  a linearized 
dynamic  model  based  on  static  wind  tunnel  data  resulted  in  a poor  pre- 
diction for  the  small-end-forward  case  and  a somewhat  better  prediction 
for  the  broadside  case.  Since  the  broadside  case  is  not  likely  to  be 
the  candidate  position  for  high-speed  operations  (due  to  drag),  a better 
analytical  method  is  needed  before  an  understanding  of  sling  load  behavior 
can  be  achieved.  Certainly,  what  has  been  learned  to  date  is  not  sufficient 
to  justify  a formal  study  of  all  the  variables  and  all  the 
possible  stabilization  schemes. 

With  this  in  mind,  our  effort  was  redirected  so  that  a dynamic  wind 
tunnel  test  could  be  performed.  The  belief  at  this  point  was  that  unsteady 
aerodynamic  effects  were  present  which  dominated  the  dynamics,  particularly 
for  the  more  important  small -end-forward  ca  e. 

The  0.1-scale  model  of  the  8-by-8-by- 20-foot  container  used  for  the  static 
tests  was  available  and  was  modified  slightly  so  that  it  could  be  hung 
in  the  tunnel  test  section.  The  suspension  system  consisted  of  two  cables 
with  end  fittings  that  allowed  the  creation  of  a number  of  single-  and 
two-point  suspension  systems.  Attach  points  were  available  at  1,  3,  5, 

7,  and  9 inches  on  either  side  of  center. 

Tufts  were  taped  to  the  model  for  airflow  identification, and  the  floor  of 
the  tunnel  was  striped  for  a camera  reference.  All  the  pertinent  variab’es 
were  measured,  including  yawing  moment  of  inertia,  weight  and  center-of- 
gravity  location,  cable  length,  and  cable  separation.  The  primary 
instrumentation  consisted  of  a 16mm  camera  located  above  the  tunnel  as 
depicted  in  Figure  1.  An  airspeed  cone  was  hung  in  the  tunnel  for  a visual 
indication  of  airspeed. 

A rigorous  analysis  of  the  model  was  made  before  the  tests  began,  using 
the  linearized  equations  of  motion  and  the  sta> ic  wind  tunnel  test  results. 
To  measure  the  yaw  damping  function  , a variety  of  cases  sensitive  to 

r 

this  parameter  were  isolated  and  incorporated  in  the  test  plan. 

The  basic  test  procedure  was  to  restrain  the  model  by  a hand-held  fork, 
bring  the  tunnel  up  to  speed,  start  the  camera,  and  then  release  the  model. 
The  initial  sway  angle  was  set  to  either  zero,  six, or  twelve  degrees. 

Zero  was  generally  used  for  the  more  active  conditions, and  six  degrees 
was  generally  used  for  the  milder  ones. 
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FIGURE  I.  0.1-SCALE  MILVAN  DYNAMIC  WIND  TUNNEL  TEST  SETUP, 


Documented  runs  totalled  47,  Including  3 calibration  runs,  17  basic 
container  runs,  22  with  fixed  fits  added,  3 with  drogue  chutes,  and  2 for 
smoke  flow  visualization  tests.  Included  in  the  basic  and  finned  tests 
were  level  and  drooped  runs  (to  vary  angle  of  attack),  narrow  cable  runs, 
unequally  spaced  cable  runs,  and  single  cable  runs. 

TEST  RESULTS 

1.  The  first  free  model  run  at  20  knots  tunnel  speed  for  the  level  box 
showed  a surprising  disparity  from  the  analog  computer  runs  using  the 
static  aerodynamic  data.  It  was  clear  that  an  energy-producing  device 
was  causing  a large  yaw  limit  cycle  of  approximately  ±20  degrees  amplitude 
at  the  bifilar  frequency.  As  the  tunnel  speed  was  increased  to  40  knots, 
the  amplitude  increased  and  the  fixed  point  of  rotation  moved  forward  to 
the  front  of  the  model.  All  stable  limit  cycle  motions  occurred  at  the 
bifilar  frequency. 

When  hung  from  a single  cable,  the  model  tur  .ed  broadside  and  developed 
a lateral  swaying  motion.  With  closely  spac  d cables,  the  center  of 
rotation  moved  toward  the  rear  of  the  model.  For  angles  of  attack  of 
-10  and  -14  degrees,  the  model  was  almost  completely  still. 

2.  Eight  finned  configurations  were  tested,  three  of  which  produced  very 
little  motion.  The  successful  configurations  incorporated  a long  fin 
projecting  out  of  the  center  of  the  rear  face  (Runs  12-15)  or  two  fins 
projecting  back  from  the  rear  comers,  simulating  opening  the  rear  doors 
of  the  MILVAN,  The  fin  length  was  20  percent  of  the  length  of  the  MILVAN. 
Fins  of  10  percent  MILVAN  length  were  not  effective.  Tilting  the  finned 
MILVAN  10  degrees  down  was  even  better,  as  the  tunnel  speed  could  be  taken 
to  67  knots.  Initial  sway  angles  of  12  degrees  were  successfully  tested 
at  this  speed. 

3.  Drogue  chutes  had  little  effect  on  the  stability  of  the  MILVAN.  One 
drogue  chute  was  simulated  by  a paint  filter  cone  and  was  attached  to  the 
center  of  the  rear  of  the  MILVAN.  This  configuration  showt 1 almost  no 
reduction  in  yaw  amplitude.  When  the  chute  was  connecteo  by  two  cables 
to  the  outer  comers  of  the  container,  the  yaw  amplitude  was  reduced  when 
no  initial  sway  angle  was  introduced;  however,  it  reached  a stable  limit 
cycle  in  yaw  after  a disturbance. 

4.  For  smoke  tests  (Runs  34  and  35),  a smoke  generator  was  hand-held 

in  the  vicinity  of  the  model  and  the  path  of  the  smoke  was  observed  from 
above  and  at  the  sides  of  ‘'he  test  section.  Smoke  released  near  the  center 
of  the  forward  face  would  billow  out  around  the  sides,  revealing  a 
separation  pattern.  As  the  model  oscillated,  the  pattern  alternated 
asymmetrically.  The  face  on  the  downstream  side  exhibited  a large  sep- 
aration plume,  while  the  plume  on  the  upstream  side  appeared  to  subside. 
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Of  more  significance  was  the  rapidity  with  which  the  plume  changed.  Smoke 
emitted  just  above  the  rear  face,  followed  one  of  two  curved  paths,  and 
appeared  to  snap  from  one  to  the  other  just  after  the  sideslip  angle  passed 
through  zero. 

5.  The  behavior  of  tufts  on  the  rear  face  was  of  considerable  interest. 
They  would  snap  back  and  forth  just  after  the  sideslip  angle  passed 
through  zero.  For  example,  as  the  rear  face  moved  to  the  left  from  a 
maximum  right  position  (nose  left),  the  tufts  would  point  left.  As  the 
rear  face  passed  through  the  neutral  position  (about  -3  degrees  sideslip), 
the  tufts  reversed  quickly  and  pointed  right.  This  action  then  reversed 
itself. 
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NONLINEAR  ANALYSIS  OF  THE  8-BY-8-BY-2Q-FQOT  CONTAINER 


The  poor  correlation  between  flight  data  and  the  results  of  the  linearized 
analysis  prompted  the  initiation  of  a nonlinear  analysis.  It  was  decided 
to  use  a small  analog  computer  so  that  on-line  comparisons  with  film 
records  could  be  made  in  real  time.  Incorporation  of  the  nonlinear 
aerodynamic  yawing  moment  and  the  extension  of  the  yaw  degree  of  freedom 
to  allow  continuous  rotation  did  not  produce  results  that  compared  with 
flight  data,  particularly  for  the  small-end-forward  case. 

INTERPRETATION  OF  DYNAMIC  WIND  TUNNEL  TEST  RESULTS 

It  was  obvious  from  the  beginning  of  these  tests  that  a dynamic  yawing 
moment  existed  that  could  not  be  attributed  to  damping  moments  due  to 
yaw  rate.  The  identification  of  this  additional  dynamic  moment  began  with 
a detailed  examination  of  the  tuft  behavior  on  the  rear  face.  It  is  likely 
that  the  "switching"  back  and  forth  of  the  tufts  was  caused  by  a rapidly 
reversing  flow  pattern  which,  in  turn,  was  caused  by  the  transition  of 
the  flow  back  and  forth  from  a partially  separated  state  to  a fully 
separated  one.  The  resulting  flow  required  time  to  establish  itself  and, 
in  the  process  of  doing  so,  generated  unsteady  flow  effects.  For  example, 
when  a two-dimensional  airfoil  is  oscillated  in  a wind  tunnel,  a vortex 
is  generated  on  the  suction  side  while  the  angle  of  attack  is  increasing. 
This  vortex  precedes  the  separation  and  is  accompanied  by  a low-pressure 
region  which  propagates  downstream  and  creates  a "pulse"  of  pitching 
moment. 

It  appears  that  this  mechanism  is  responsible  for  creating  the  dynamic 
yawing  moment.  Specifically,  as  sideslip  reaches  about  3 degrees  and  is 
increasing,  a vortex  propagates  downstream  along  the  downstream  side. 

When  it  reaches  the  rear  face  where  the  vortex  can  now  "draw  in"  more  flow, 
it  does  so,  and  a small  sharp  pulse  of  negative  pressure  appears  on  that 
side  of  the  rear  face.  This  causes  the  tufts  to  reverse  and  creates  a 
small  pulse  of  yawing  moment  that  reinforces  the  oscillation.  With  this 
hypothesis  in  mind,  Figure  2 depicts  the  dynamic  yawing  moment  evolved 
from  the  analog  matching  study.  Note  that  a moment  hysteresis  is  present; 
i.e.,  a dynamic  moment  is  present  only  when  the  product  of  sideslip  and 
sideslip  rate  is  positive  and  sideslip  remains  within  certain  limits. 

ANALOG  MATCHING  STUDY 


The  motion  of  the  tufts  and  the  smoke  pattern  suggested  that  a vortex 
was  being  shed  as  the  model  rotated  through  nearly  zero  sideslip.  The 
purpose  of  the  analog  matching  study  was  to  find  a vortex  model  that  could 
be  added  to  the  wind  tunnel  static  aerodynamics  and  conventional  rotary 
derivatives  which  Would  satisfactorily  match  *-he  observed  dynamic  tunnel 
test  results  over  a wide  range  of  conditions.  Ground  rules  were  established 
for  the  vortex  model  as  follows: 
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MODELED  FROM  THE  RESULTS 
OF  A DYNAMIC  WIND  TUNNEL  TEST 
(STROUHAL  NUMBER  * 0.1  RADIAN. 
BASED  ON  WIDTH  b,  ) 


FIGURE  2.  DYNAMIC  (UNSTEADY)  YAWING  MOMENT  FUNCTION. 


1.  Tue  model  should  be  compatible  with  established  concepts  of  unsteady 
flow  effects. 

2.  No  changes  to  accommodate  a match  at  different  speeds  should  be 
necessary. 

3..  Changes  could  be  made  to  accommodate  different  angles  of  attack  and 
the  addition  of  fins. 

The  matching  conditions  consisted  of  first  generating  an  oscilloscope 
displ  y representing  the  centerline  of  the  MILVAN,  scaled  so  that  the 
sway  and  yaw  modes  appeared  as  if  viewed  from  above.  Next  to  the 
oscilloscope,  the  films  of  the  tests  were  projected  at  approximately  the 
same  size.  A compute  * , programmed  with  the  appropriate  initial  conditions, 
was  set  to  start  computing  when  the  projected  film  showed  the  model  being 
released. 

It  was  obvious  at  the  start  that  an  exact  movement-for-movement  match 
could  not  be  made;  therefore,  an  attempt  was  made  to  match  the  limit  cycle 
frequence,  growth  and  decay  rates,  and  fixed  point  of  rotation  or  phasing 
of  the  modes.  Since  was  an  unknown  parameter,  initial  attempts  were 

r 

made  to  match,  using  this  parameter.  It  was  found  that  a match  could  be 

made  at  various  levels  of  C_,  . Hence,  it  was  set  to  a "reasonable"  level 

N 

r 

and  not  changed.  The  larger  the  CM  , the  larger  the  value  of  the  unsteady 

r 

moment  that  was  required  to  match  the  wind  tunnel  data.  This  led  to  the 
energy  balance  concept  presented  in  a following  subsection.  The  shape 
and  magnitude  of  the  unsteady  moment  function  had  considerable  influence 
on  the  growth  and  steady-state  magnitude  of  the  limit  cycle.  The  shape 
could  not  be  defined  precisely  because  of  the  inaccuracy  of  the  test 
data  and  the  use  of  an  eyeball  match;  however,  as  more  test  runs  were 
matched,  more  confidence  was  gained  in  the  model. 

Some  of  the  conditions  that  have  a substantial  impact  on  the  shape  of 
the  unsteady  function  (Figure  2)  are  illustrated  by  the  following  examples. 
Three  runs  at  a 14-degree  angle  of  attack  displayed  an  extreme  speed 
sensitivity  and  a dual  mode  stability.  These  cases  and  others  tended 
to  have  a small-amplitude  stable  limit  cycle  with  no  initial  disturbance 
and,  when  disturbed,  a large-amplitude  limit  cycle.  An  increase  in 
speed  in  one  run  at  -14  degrees  of  attack  angle  caused  a change  in  limit 
cycle.  There  were  no  cases  of  limit  cycle  amplitudes  between  approximatel’' 
±2  and  ±10  degrees.  These  conditions  were  met  by  reducing  the  magnitude 
of  the  unsteady  function  at  sideslip  angles  near  zero  and  forcing  it  to 
zero  beyond  ±9.5  degrees  sideslip.  After  this  form  of  unsteady  moment 
function  was  used,  only  changes  in  magnitude  were  required  to  match  the 
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various  configurations.  The  fact  that  the  unsteady  function  was  zero 
beyond  ±10  degrees  sideslip  and  the  fact  that  the  container  had  its  best 
stability  at  angles  of  attack  of  ±10  degrees  are  probably  related. 

Because  of  the  phasing  of  the  sway  and  yaw  motions,  one  point  on  the 
container  appeared  to  be  fixed  relative  to  the  tunnel.  This  is  referred 
to  as  the  "fixed  point  of  rotation."  It  is  interesting  to  note  that  this 
fixed  point  of  rotation  was  near  the  front  and  moved  forward  with  increased 
speed  when  X = 9 inches.  However,  with  reduced  cable  separation 
(X^  = 2 inches),  the  fixed  point  was  near  the  rear  of  the  container. 


In  all  cases  with  the  wide  cable  separation,  the  sway  mode  quickly  dis- 
sipated and  the  sway  motion  was  forced  by  the  yaw  mode.  With  the  cable 
separation  such  that  the  yaw  frequency  and  sway  frequency  are  nearly  the 
same,  the  modes  couple  and  instability  of  the  classical  nature  results. 

The  single-cable  cases  matched  to  near  perfection. 

The  analog  model  matched  quite  well  with  the  wind  tunnel  results.  The 
limit  cycle  frequency  matched  to  approximately  1/4-cycle  in  10  cycles, 
or  2.5  percent.  Nothing  was  done  to  the  model  to  change  the  frequency 
because  it  was  primarily  controlled  by  the  cable  separation.  With  small- 
amplitude  limit  cycles  the  static  yawing  moment  had  an  effect;  however, 
with  large-amplitude  limit  cycles  it  averaged  out  to  nearly  zero.  The 
growth  rate  of  the  limit  cycle  is  probably  accurate  to  between  10  and  20 
percent.  The  amplitude  of  the  limit  cycle,  while  the  most  easily  changed 
variable,  is  probably  15  to  20  percent  in  error,  due  primarily  to  the  camera 
angle  not  being  vertical. 

The  shape  of  the  unsteady  moment  curve  is  undoubtedly  in  gross  error 
because  of  the  lack  of  precise  data.  The  straight  lines  of  the  function 
do  not  reflect  the  lags  associated  with  unsteady  aerodynamics;  however, 
the  general  characteristics  of  the  unsteady  moment  which  evolved  must  have 
some  correlation  with  reality  because  the  predicted  MILVAN  motions  agreed 
well  over  a wide  range  of  speeds  and  geometric  conditions. 

DYNAMIC  MODEL  AND  PREDICTIONS 

Figure  3 is  a block  diagram  of  the  dynamic  model.  A logic  representation 
of  the  unsteady  yawing  moment  function  is  shown  and  the  static  yawing 
moment  is  also  included.  The  model  has  a single  degree  of  freedom  describing 
the  motions  about  the  yaw  axis  only,  which  is  usually  a stable  limit  cycle. 
This  model  is  valid  as  long  as  the  yawing  limit  cycle  frequency  (bifilar 
frequency)  is  significantly  greater  than  the  sway  pendulum  frequency. 

If,  for  example,  the  bifilar  frequency  is  lowered  to  about  1.5  times  the 
sway  pendulum  frequency,  the  sway  mode  will  be  unduly  perturbed  by  yaw 
and  a classical  instability  usually  begins  to  appear. 
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FIGURE  3.  BLOCK  DIAGRAM  OF  THE  DYNAMIC  MODEL  OF  THE  8-BY-8-BY-20-FOOT 
CONTAINER  FOR  THE  LIMIT  CYCLE  CLASS  OF  MOTION. 


ENERGY  BALANCE 

PREDICTION  OF  LIMIT  CYCLE  (DYNAMIC  MODEL) 

An  attempt  was  made  to  estimate  the  amplitude  of  the  yaw  limit  cycle  based 
on  the  conservation  of  energy  concept.  Specifically,  the  yaw  motion  is 
assumed  to  be  a stable  limit  cycle  at  the  bifilar  frequency,  and  the 
energy  added  to  the  system  by  the  unsteady  aerodynamics  is  exactly  removed 
by  the  damping  function  C^  with  the  assumption  of  sinusoidal  motion. 

r 

Integration  of  the  unsteady  moment  function  into  energy  is  complex  for 
amplitudes  less  than  ±?.5  degrees;  however,  for  larger  amplitudes,  a 
closed-form  solution  can  be  derived  as  given  below. 

= (|Y  I > .166  rad) 

max  ttCn  bm  '‘maxi  “ 

r * 

In  this  expression,  the  negative  sign  reflects  the  damping  offered  by 
CN  , which  is  usually  taken  as  minus;  0.1126  is  the  integral  value  of  the 
r 

unsteady  moment  function,  and  C is  the  normalized  value  of  the  unsteady 
function  and  is  a function  of  fins  and  angle  of  attack.  For  a peak  amplitude 
greater  than  0.166  radian,  C is  only  a function  of  angle  of  attack  as 
given  in  Figure  4.  A straight-line  function  was  assumed  between  the 
measured  points. 


FIGURE  4.  EFFECT  OF  ANGLE  OF  ATTACK  AND  FINS. 
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It  is  convenient  to  express  the  relation  in  terms  of  the  Strouhal 

number  (S„  = mb/V  ) thus: 

1 7 ° 


_ /-2&.&.1226  nr 
V "S.  V st 

where  b is  container  width,  <«  is  the  bifilar  frequency  in  rad/sec,  and 

7 

V is  forward  speed  in  ft/sec.  The  value  of  C.,  derived  from  the  analog 

O 

matching  study  is  -0.94.  When  this  is  substituted  in  the  equation,  a 
relatively  simple  expression  results: 


- VSS*  VV  °-2f 


Figure  5 shows  the  model  data  points,  the  full-scale  points  from  Reference 
2,  computer  predictions,  and  the  energy  balance  estimate.  The  agreement 
between  the  computer  results  and  the  energy  balance  prediction  is  very 
good.  The  flight  data  from  Reference  2 also  shows  good  agreement  with  the 
model  for  short  cables  and  low  yaw  amplitude  (±15  degrees).  The  agreement 
between  the  energy  balance  prediction  and  flight  data  for  large- 
amplitude  (long  cables)  yawing  is  not  as  good.  The  reason  for  this  is 
that  the  simple  formula  assumes  the  motion  to  be  sinusoidal  and  the 
computer  prediction  shows  it  to  be  approaching  a triangular  form.  This 
will  certainly  affect  the  calculation  for  the  damping  function  CN  . The 

r 

reasonably  good  match  between  the  model  prediction  of  flight  results  and 
the  flight  results  themselves  is  remarkable  considering  the  strong  effect 
of  angle  of  attack  and  the  lack  of  documentation  of  angle  of  attack  in 
the  flight  data. 

STROUHAL  NUMBER 

This  is  a nondimens ional  frequency  used  in  dynamic  wind  tunnel  testing. 

It  is  a frequency  in  radians  per  second  or  cycles  per  second  multiplied 
by  length  and  divided  by  velocity.  Holding  this  number  constant  from 
model  to  full  scale  gives  the  same  result  as  for  Froude's  number;  namely, 
full-scale  speed  is  3.16  times  model  speed.  The  Strouhal  number  is  useful 
in  describing  vortex  shedding  and  other  properties  of  a turbulent  wake. 

It  is  reasonable  to  hold  it  constant  for  dynamic  wind  tunnel  testing  of 
this  kind  because  in  doing  so  the  relation  between  wake  motions  (vorticity, 
etc.)  and  model  motions  is  maintained,  thus  increasing  the  probability  of 
correctly  representing  the  unsteady  aerodynamic  effects. 
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FIGURE  5.  COMPARISON  OF  COMPUTER  MODEL  1 ATA  WITH  FLIGHT  RESULTS 


MODEL  347  HELICOPTER  SIMULATION 


SIMULATION  FACILITY 


The  Northrop  simulation  facility  used  for  the  Model  347  helicopter  simu- 
lation comprises  a large-amplitude  beam-type  motion  subsystem  that  carries 
a cockpit,  projectors  for  sky/earth/horizon/target  imaging,  and  a spherical 
display  screen  for  the  visual  display,  all  on  the  end  of  a 30-foot  beam 
as  shown  in  Figure  6.  This  motion  system  provides  vertical  and  lateral 
beam  translation  and  cockpit  pitch,  roll,  and  yaw  rotation  as  shown  in 
Table  4. 


TABLE  4.  SIMULATOR  PERFORMANCE 


COCKPIT  MOTION 
AMPLITUDES 

DISPLACEMENT 
(USABLE  TRAVEL) 

PEAK 

NO-LOAD 

VELOCITY 

PEAK 

ACCELERATION 

Vertical  Translation  (Beam) 
Lateral  Translation  (Beam) 
Pitch  Rotation  (Cockpit) 
Roll  Rotation  (Cockpit) 

Yaw  Rotation  (Cockpit) 

± 10.4  Ft* 
± 10.1  Ft* 
± 25  Deg 
± 20  Deg 
± 20  Deg 

13.2  Ft/Sec 

9.0  Ft/Sec 

1.0  Rad/Sec 

1.1  Rad/Sec 
1.3  Rad/Sec 

2.5  G 

1.6  G 

16.0  Rad/Sec _ 

4.6  Rad/Sec 

8.6  Rad/Sec 

■-Measured  at  the  pilot's  eye  level 

In  operation,  the  motion  subsystem  produces  movements  in  precise  phase  and 
amplitude,  with  electrical  signals  supplied  by  an  integrated  hybrid 
computer  complex  to  provide  the  pilot  with  realistic  perceptual  cues  re- 
lated to  actual  motions  of  the  simulated  aircraft.  The  large  vertical  and 

[lateral  translations  of  the  30-foot  gimbaled  beam  are  controlled  by  two 

electrohydraulic  servo  systems.  The  cockpit  is  also  gimbaled  on  the  beam 
to  provide  pitch,  roll,  and  yaw  rotations,  all  controlled  by  electro- 
hydraulic  servos.  Figure  7 shows  the  gimbaled  structure  supporting  the 
| cockpit.  The  dynamic  performance  of  the  large-amplitude  simulator  is 

| shown  graphically  in  Figure  8. 
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AND  BEAM  END  CONFIODRA  I' (ON 


FIGURE  8.  LARGE -AMPLITUDE  SIMULATOR  - DYNAMIC  PERFORMANCE 


The  visual  display  subsystem  consists  of  two  different  projectors  mounted 
above  the  cockpit  inside  the  spherical  housing.  A sky-earth  projector 
provides  the  pilot  with  a horizon  reference  that  continually  gives  visual 
indication  of  the  simulated  aircraft  attitude.  A target  projector  is  used 
primarily  for  air-to-air  simulation  but  can  also  be  used  for  air-to-ground 
work.  This  target  information  is  projected  into  the  screen  from  a cathode- 
ray  tube  located  above  and  behind  the  pilot's  head. 

The  force  feel  for  primary  controls  is  provided  by  a three-axis  electro- 
hydraulic  servo  system.  This  system  consists  of  a conventional  two-axis 
center  stick  and  rudde . pedals.  Force  is  applied  to  the  controls  by 
rotary  hydraulic  actuators,  and  the  force-feel  characteristic  is  provided 
by  the  console  servo  electronics  and/or  the  computer.  A unique  capability 
of  this  force  servo  system  is  that  it  is  capable  of  stable  operation  with 
zero  force  applied  to  the  controls,  i.e.,a  true  feel  stick. 

The  monitor  and  control  console  is  the  central  operating  station  for 
simulator  control,  monitoring,  test,  and  maintenance.  The  console  is  shewn 
in  the  lower  part  of  Figure  9.  It  houses  the  electronic  portions  of  the 
servo,  cockpit,  and  television  equipment  and  provides  the  signal  interface 
between  the  simulator  and  the  computer.  An  intercom  system  in  the  console 
provides  communications  between  console  operator,  pilot,  and  the  computer 
area. 

Other  support  components  of  the  simulator  facility  include  an  aircraft- 
engine  sound  generator,  an  anti-G  suit  pressurization  system,  a 3000-psi 
hydraulic  power  source,  a hybrid  computer  complex,  and  software  programs 
that  drive  the  simulator. 

MODEL  347  HELICOPTER  SIMULATION  SETUP 


To  understand  the  effects  of  various  sling  load  stabilization  concepts, 
a sophisticated  simulation  of  the  Model  347  helicopter  was  developed  and 
mechanized  on  Northrop' s advanced  large-amplitude  simulator.  Linearized 
equations  of  motion  for  the  sling  load  were  used, and  the  dynamic  model  of 
the  8-by-8-by-20-foot  cargo  container  (MILVAN)  was  incorporated  in  the 
simulation. 

The  existing  simulator  c >ckpit  was  modified  to  resemble  the  right-hand 
seat  of  the  Model  347  helicopter.  An  instrument  panel  was  fabricated  to 
duplicate  the  configuration  of  the  Model  347  helicopter;  all  primary 
instrumentation  was  actively  driven  by  the  computer.  Instrumentation 
included  attitude  director  indicator,  horizontal  situation  indicator, 
airspeed,  radar  and  barometric  altimeters,  rate  of  climb,  engine  rpm 
and  rotor  torque  indicators.  A Boeing -supplied  collective  control  and 
Automatic  Flight  Control  System  (AFCS)  control  panel  were  also  installed 
in  the  simulator  cockpit  and  mechanized  to  operate  similar  to  that  in 
the  helicopter.  A computer  mechanization  of  the  control  force-feel  system 
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was  developed  to  provide  control  forces  as  a function  of  airspeed,  pitch 
rate,  AFCS  drive  commands,  and  the  magnetic  brake  function. 

Additional  instruments  added  to  the  cockpit  for  performance  of  the  hover 
tasks  were  longitudinal  and  lateral  airspeed  indicators  and  a computer- 
generated CRT  display.  The  CRT  display  could  be  presented  in  three  modes: 

1)  Aircraft  motion  relative  to  a fixed  point  on  the  ground. 

2)  Motion  of  a fixed  point  on  the  ground  relative  to  the  aircraft. 

3)  Slung  load  motion  relative  to  the  aircraft. 

The  scale  of  the  CRT  display  could  be  changed  by  the  pilot  to  depict 
relative  position  (a  coarse  scale  of  500  ’et  to  a fine  scale  of  5 feet 
in  position). 

A validation  effort  to  establish  confidence  in  the  simulation  preceded  the 
principal  test  phase.  This  validation  consisted  of  setting  in  the  coeffi- 
cients which  produced  a simulation  of  a known  flight  condition  and  then 
having  the  pilot  who  flew  the  knovai  flight  condition  fly  the  simulator. 
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VALIDATION 


The  mathematical  model  validation  of  the  347  Advanced  Tandem-Rotor 
Helicopter  was  accomplished  by  comparison  of  steady-state  trim  and 
dynamic  response  data  with  theoretical  predictions. 

Level- flight  data  were  obtained  at  airspeeds  ranging  from  40  knots 
rearward  to  180  knots  forward  at  a gross  weight  of  45,000  pounds  and  a 
center  of  gravity  (c.g.)  condition  of  12  inches  aft  of  the  nominal  c.g. 
location. 

Trim  validations  were  performed  using  two  different  approaches  to  generate 
rotor  force  and  moment  parameters:  the  "classical"  equation  approach  and 
the  "rotor  map"  approach.  The  classical  approach  utilized  an  equation 
format  to  generate  rotor  forces  and  moments,  whereas  the  rotor  mapping 
technique  involved  the  use  of  precalculated  values  for  thrust,  propulsive 
force,  and  power.  These  precalculated  values  were  then  stored  in  the 
PDP9B  digital  computer  which  was  used  as  a function  generation  device 
for  the  math  model  simulation. 

The  correlation  between  values  predicted  by  Boeing  Vertol  (B-29)  and 
Northrop  hybrid  results  of  steady-state  trim  data  is  shown  on  Figure  10. 

Data  agreement  as  shown  on  these  figures  is  generally  good  for  all  trim 
conditions.  Improvements  due  to  the  "rotor  map"  approach  to  rotor  forces 
and  moments  are  noticeable  for  the  higher  airspeeds  (140  to  180  knots). 

Time  history  responses  were  obtained  for  both  the  augmented  and  unaug- 
mented flight.  The  augmented  flight  consisted  of  longitudinal,  lateral, 
and  directional  stability  augmentation  systems.  Control  inputs 
consisting  of  a series  of  steps  and  pulses  were  introduced  into  the 
longitudinal,  lateral,  and  directional  axes.  The  resulting  dynamic 
characteristic  comparisons  with  expected  results  were  good. 

The  validation  of  the  sling  load  model  was  accomplished  with  a 15,000- 
pound  external  load  and  a 30,000-pound  helicopter  at  12  inches  aft  c.g. 
at  standard  sea  level  conditions.  Sling  load  pendulum  (longitudinal 
anu  lateral)  and  bifilar  (yaw)  modes  were  excited  by  separately  inputting 
1-inch  and  ^-second  pulses  into  the  longitudinal,  lateral,  and  directional 
axes  of  the  347  analog  control  system.  Cable  angle  responses  were  measured 
and  frequency  and  damping  information  recorded.  The  results  obtained 
correlated  «ell  with  the  predictions. 

The  time  history  results  of  the  aircraft  and  sling  load  models  as  well 
as  the  static  trim  data  were  compared  with  similar  data  acquired  from  the 
Boeing  Vertol  simulation  of  the  347  aircraft.  The  availability  of  data 
from  the  Boeing  Vertol  simulation  enhanced  the  validation  of  the  mathematical 
model  at  Northrop.  Data  correlations  from  the  two  math  models  were 
exact. 
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Airspeed  - knots 


Figure  10.  NORTHROP  HYBRID  VALIDATION  OF  THE 

MODEL  347  HELICOPTER  (SHEET  1 OF  10). 

39 


Degree* 


I 


a 

a. 

O 

<i> 


OC 

O 

Oi 


Gross  weight  * 46.000  lb 
CG  • 12  in.  aft 
S<a  level 


VT  * 728  ft/sec 


B 29  - 5366.8186. 8373 
Northrop  hybrid 
O Classical 
A Lmt  maps 


Airspeed  - knots 


Figure  10.  (Sheet  5 of  10) 
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The  Boeing  Vertol  math  model  was  validate-!  by  comparing  theoretical 
results  and  model  simulation  results  with  actual  flight  test  data. 

The  comparison  of  the  Vertol  hybrid  model  with  flight  test  data  showed 
good  correlation  in  the  speed  regimes  modeled.  The  results  of  the  Vertol 
simulation  of  the  347  aircraft  were  then  used  as  the  basis  for  comparison 
in  the  validation  of  the  347  aircraft  math  model  here  at  Northrop. 

Upon  completion  of  the  mathematical  model  validation  the  LAS/WAVS  simulator 
validation  was  performed,  utilizing  a subjective  approach  in  which  pilot 
opinion  was  the  sole  indicator  of  fidelity.  With  a test  pilot  from  the 
Boeing  Vertol  347  Aircraft  Flight  Test  Group  flying  the  simulator,  changes 
in  the  interface  coefficients  were  made  to  "tune"  the  simulator  such 
that  the  motion  and  visual  stimuli  presented  to  the  pilot  resembled  those 
experienced  on  the  actual  347  aircraft. 
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SIMULATION  RESULTS 


MASS  RATIO  EFFECTS 


A two-degree-of-f reedom  nonlinear  analysis,  while  rigorous,  does  not 
predict  the  limit  cycle  amplitude  with  precision  when  the  load  is  attached 
to  a helicopter  of  finite  mass.  For  the  empty  MILVAN  with  spreader,  the 
ratio  of  load  mass  to  total  mass  is  8705  (37,000  + 8705)  = 0.19.  The  ratio 
of  load  yawing  inertia  to  helicopter  inertia  is  even  smaller,  9500/367,000  = 
0.026.  It  would  seem  that  the  condition  of  infinite  helicopter  mass  would  be 
satisfied  for  the  purpose  of  employing  the  two-degree-of-f reedom  analysis. 
However,  when  the  dynamic  model  was  incorporated  into  the  simulation, the 
limit  cycle  amplitude  matched  the  earlier  model  results  precisely, 
but  only  if  the  helicopter  was  held  fixed  (by  manipulating  the 
computers. ) 

When  the  helicopter  was  released  after  trimming,  the  limit  cycle  amplitude 
was  reduced  to  0.75  of  the  helicopter  fixed  value.  The  difference  amounted 
to  about  three  degrees.  The  initial  flight  data  for  the  validation  flight 
condition  was  listed  as  having  an  amplitude  of  ±15  degrees,  but  the  axes 
reference  was  not  stated.  In  the  simulation,  the  reference  is  body  axes, 
since  the  helicopter  yaws  at  the  limit  cycle  frequency,  the  load  yaws 
±12  degrees  relative  to  the  helicopter,  and  the  helicopter  yaws  ±3  degrees 
relative  to  the  earth.  If  the  phasing  between  them  is  180  degrees,  the 
load  yaw  angle  relative  to  the  helicopter  is  ±15  degrees.  This,  then, 
identifies  a minor  inconsistency  in  the  comparison  of  the  model  results 
and  flight  simulation  with  flight  data. 

Figure  11  shows  the  results  of  an  attempt  to  predict  the  carrying  speed 
boundary  for  the  empty  container  and  spreader  beam,  using  the  model  and  the 
simulation-determined  limit  cycle  restrictions.  The  yaw  angle  is 
referenced  to  helicopter  body  axes.  The  0.75  factor  was  applied  to  the 
model  prediction  of  yaw  amplitude.  The  model  prediction  and  the  carrying 
speed  amplitude  restriction  curves  intersect  at  60  knots.  The  flight  report 
indicates  that  the  pilots  would  not  recommend  carrying  this  load  beyond 
65  knots.  This  is  a very  good  agreement,  and  it  can  be  concluded  that  the 
dynamic  model  is  a good  representation  of  the  real-world  effects, 
particularly  when  integrated  into  a sophisticated  simulation.  However, 
use  of  the  dynamic  model  alone  to  predict  carrying  speed  boundaries  is 
subject  to  the  0.75  uncertainty;  i.e.,  the  effect  of  load-to-helicopter 
mass  ratio. 

EMPIRICAL  DATA  FIT 

When  confidence  in  the  simulation  was  established,  a series  of  tests  wis 
conducted  to  determine  the  limit  cycle  amplitude  that  the  pilots  felt  was 
the  maximum  allowable.  To  conduct  this  test  series,  the  cable  separation 
was  varied  so  that  the  bifilar  frequency  could  be  controlled.  Five  data 
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points  were  obtained  lot  the  tangc  ot  oiriLr  iicqutnv.  tim<i  1.  Lad/sec 
to  5.22  rad/sec.  For  the  spreader  beam  configuration  with  7.5-foot  cables, 
this  corresponds  to  a cable  separation  range  from  7.  7 feet  to  ?•)  feet, 
respectively. 

The  five  data  points  revealed  that  for  a given  airspeed  the  pilot  will 
accept  large -amplitude  limit  cycle  load  otions, provided  these  motions  are 
slow;  i.e.,  low  bifllar  frequency.  Before  he  will  carry  the  load  any  faster, 
the  load  motion  amplitude  must  be  made  lower  or  its  frequency  reduced. 


Curiously  enough,  these  five  data  points  can  be  fitted  very  well  by  the 
simple  relation 


_ -091 


max 


This  form  resulted  from  an  estimation  of  limit  cycle  amplitude,  considering 
that  the  energy  added  to  the  system  by  the  unsteady  aerodynamics  was  equal 
to  that  removed  by  the  damping  function. 

This  relation  is  clarified  by  the  plot  of  Figure  12,  in  which  the  maximum 
allowable  amplitude  is  plotted  versus  airspeed  and  bifilar  frequency. 

The  five  simulation  data  points  are  also  shown. 


AMPLITUDE  AND  FREQUENCY. 


STABILIZATION.  CONCEPTS 


It  has  been  demonstrated  that  use  of  the  dynamic  model  with  simulation 
gives  credible  results.  That  is,  an  accurate  assessment  of  carrying  speed 
limitation  can  be  determined  for  all  combinations  of  cable  geometry  and 
load  mass. 

We  have  also  seen  how  a carrying  speed  limitation  based  on  pilot  judgment 
was  evolved.  The  limitation  has  some  relation  to  a simple  energy  balance 
expression  which  describes  the  load  limit  cycle  amplitude.  The  underlying 
laws  governing  the  limitation  criteria  for  the  limit  cycle  are  not  well 
understood  at  this  time.  While  it  is  important  to  acquire  this  understanding, 
it  was  more  important  for  the  purpose  of  this  study  to  determine  what  factors 
affect  the  carrying  speed  limitation  and  how  it  can  be  extended  within 
the  current  state  of  the  art.  To  this  end,  the  following  stabilization 
concepts  to  extend  carrying  speed  were  examined: 

1.  Optimizing  the  suspension  system. 

2.  Maneuvering  the  helicopter  to  take  advantage  of  the  inherently 
good  aerodynamic  properties  of  the  bare  container. 

3.  Modifying  the  container  shape  to  improve  its  aerodynamic 
characteristics. 

4.  Using  an  active  stability  augmentation  system  that  generates  moments 
through  the  cables  by  moving  cable  guides  in  the  helicopter. 

5.  Modifying  the  helicopter  stability  augmentation  system  to  include 
load  motion  feedback. 

6.  Using  an  active  stability  augmentation  system  employing  aerodynamic 
moments  derived  from  an  airfoil  or  similar  device. 

A summary  of  key  properties  of  each  concept  is  given  in  Figure  13. 
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Optimized  suspension:  two-point, 

pendulum  length  less  than  20  ft  and 
20-ft  parallel  cable  separation. 

Aft  winch  to  optimize  angle  of 
attack  to  -10  degrees.  Sensitive 
to  piloting  technique. 
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Optimum  maneuvering;  carrying 
maneuver  is  designed  to  maintain 
an  angle  of  attack  or  sideslip  of 
10  degrees,  while  minimizing 
mission  time.  Contingency  man- 
euvers in  case  of  emergency 
(collision  avoidance,  etc) 
require  additional  pilot  training. 


WAKE  vo  ^ 
SPLITTER  ^ 


<0 


ROUNDED 
NOSE 


'0 

V 


CORNER 
DEFLECTORS3 


Modifying  container  shape;  wake 
splitter  plate,  aft  fins,  comer 
deflector,  round  nose  all  tend  to 
reduce  unsteady  aerodynamic  yaw 
moments.  Sensitive  to  normal 
operational  abuse. 
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CABLE 
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ACTUATORS 


Cable  guide  actuators;  two  cable  guide 
actuators  deflect  the  cables  laterally 
(collectively)  in  response  to  load  sway 
rate,  and  laterally  (d*ff.)  in  response  to 
load  yaw  rate  (relative  to  heli.).  Desen- 
sit'  zes  the  problem,  permits  brisk  maneu- 
vering, but  requires  additional  actuator 
and  load-motion-sensor  hardware. 
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Helicopter  SAS  inputs  from  load 
motion  sensors;  lateral  sway  rate* 
lateral  SAS;  load  yaw  rate-»yaw 
SAS.  Improves  carrying  speed 
slightly,  degrades  helicopter 
ride  qualities. 


ET 


AUXILIARY 
AIRFOIL 
DRIVEN  BY 
YAW  RATE 
SENSOR 


Auxiliary  airfoil;  driven  by  load 
yaw  rate?  May  improve  carrying 
speed  if  airfoil  operates  effec- 
tively In  wake  of  container;  air- 
foil may  possibly  be  placed  forward 
of  load.  Complex  to  handle,  requires 
additional  hardware. 


FIGURE  IX  KEY  PROPERTIES  OF  STABILIZATION  SYSTEMS  STUDIED. 
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The  effects  of  load  mass  have  significant  impact  on  sling  load  carrying 
speeds.  One  of  the  tests  performed  during  the  Model  347  helicopter  simu- 
lation was  an  evaluation  of  a 15,000-pound  container  on  a two-point  suspension 
carried  by  a 30,000-pound  347  helicopter.  The  total  mass  was  nearly  the 
same  as  for  tests  with  the  empty  container,  but  the  carrying  speed  limitation 
was  higher  for  this  configuration.  The  load  oscillated  in  yaw  as  before, 
but  the  maximum  load  yaw  angles  were  much  less,  allowing  the  evaluation 
pilot  to  reach  simulator  speeds  up  to  150  knots.*  The  worst  case  occurs 
with  the  empty  container. 

SUSPENSION  SYSTEM  OPTIMIZATION 


When  single-point  suspension  is  used,  the  effective  pendulum  length  can 
rarely  be  made  much  shorter  than  the  length  of  the  container.  Under  these 
conditions,  an  acceptable  speed  of  40  knots  is  usually  recommended.  At 
this  speed,  the  container  trails  in  a small-end-forward  position  at  a low 
angle  of  attack  in  level  flight.  To  take  advantage  of  the  container's 
angle  of  attack,  use  of  slings  of  different  length  would  permit  the  load 
to  be  rigged  in  a nose-down  attitude.  An  initial  rigging  angle  of  about 
five  degrees  may  reduce  the  rapid  nose  departure  from  the  small-end-forward 
position,  thus  allowing  the  pilots  to  increase  their  speed,  perhaps  to  50 
knots.  This  configuration,  however,  will  demonstrate  high  sensitivity  to 
climbs  and  descents,  since  load  angle  of  attack  will  vary  widely  during 
these  maneuvers. 

There  is  a greater  possibility  for  improvement  using  suspension  from  two 
or  more  points.  The  bifilar  and  pendulum  frequencies  can  be  separated  to 
eliminate  classical  sway  mode  instabilities.  For  the  empty  container, 
this  requires  a parallel  cable  separation  of  at  least  12  feet,  considering 
that  in  practical  sling  arrangements  the  bifilar  length  is  generally 

less  than  the  pendulum  length  jJ^.  Increasing  cable  separation  to  the  maximum 

practical  extent  is  also  advisable.  Furthermore,  the  container  may  be 
rigged  so  that  the  load  angle  of  attack  is  about  -10  degrees  under  level 
cruise  conditions.  For  the  container  and  spreader  beam  carried  by  a 
37,000-pound  Model  347,  a nose-down  rigging  angle  of  two  degrees  produces 


*It  is  doubtful  that  the  actual  Model  347  helicopter  could  carry  a 15,000- 
pound  container  faster  than  120  knots.  It  could  be  reached  in  the  simulator 
because  the  power  plant  limitations  were  artificially  removed  for  these 
tests. 
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this  condition  at  180  knots.  This  angle  will  vary  slightly  with  different 
helicopters.  A further  possibility  is  to  winch  the  load  up  so  that  the 
cable  length  is  zero.  The  theoretical  bifilar  frequency  is  infinite 
as  the  yaw  stiffness  is  now  provided  by  the  triangular  sling-box  arrangement. 
The  yaw  stiffness  will,  in  practice,  be  a large  but  finite  value  and  this 
should  allow  higher  speeds.  There  is  a danger,  however,  that  the  increased 
bifilar  frequency  will  now  couple  with  the  Von  Raman  vortex  shedding 
frequency.  If  this  happens,  the  problem  may  return,  but  in  a much  higher 
frequency  range,  where  vibration  loads  may  be  limiting.  It  is  also  not 
certain  that  the  slings  will  always  remain  taut  for  this  configuration. 

MANEUVERING  THE  HELICOPTER 


Maneuvering  the  helicopter  to  take  advantage  of  the  inherently  good  aero- 
dynamic properties  of  the  bare  container  is  a good  way  to  extend  allowable 
carrying  speeds.  During  the  simulation  of  the  empty  container  with  spreader, 
a curious  ambiguity  appeared.  When  the  container  was  rigged  level,  a speed 
limitation  of  90  knots  would  be  voiced  by  the  evaluation  pilot, but  on  sub- 
sequent runs  he  was  able  to  reach  a speed  of  165  knots.  Upon  closer 
examination,  it  was  found  that  on  the  higher  speed  run  he  had  accelerated 
through  the  region  of  90  knots.  Obviously,  the  nosedown  attitude  during 
acceleration  caused  the  load  to  assume  an  angle  of  attack  which  reduced 
the  limit  cycle  amplitude.  Figure  9 shows  a region  between  60 
and  90  knots  where  no  operation  is  expected.  However,  with  an  additional 
angle  of  attack,  the  two  curves  could  not  intersect  until  a much  higher 
speed  was  reached.  This  point  emphasizes  the  possibility  of  following 
a flight  profile  that  minimizes  the  limit  cycle  amplitude.  This  could 
be  achieved  by  climbing,  descending,  or  introducing  sideslip,  recalling 
that  a sideslip  of  10  degrees  or  more  suppresses  the  limit  cycle.  These 
properties  make  this  case  sensitive  to  pilot  strategy,  but  of  more  importance 
is  the  possibility  that  a series  of  maneuvers  could  be  followed  that  would 
allow  the  attainment  of  power-limited  airspeeds  for  this  container. 
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MODIFYING  THE  CONTAINER  SHAPE 


Modifying  the  container  shape  to  improve  its  aerodynamic  properties  is  a 
powerful  but  less  practical  way  of  extending  carrying  speed.  Specifically, 
it  was  found  that  the  addition  of  a wake  splitter  plate  or  rear  fins 
suppressed  the  unsteady  moments  as  much  as  the  introduction  of  the  optimum 
10-degree  angle  of  attack.  The  one  configuration  tested  with  fins  could 
be  carried  to  a speed  of  120  knots,  as  opposed  to  90  knots  without  them. 
Requirements  for  the  fins  and  the  wake  splitter  plate  are  modest.  The  fine 
could  be  the  aft  doors,  <.nd  the  splitter  could  be  a plate  about  10  feet 
in  length  mounted  on  the  centerline  of  the  aft  face.  These  configurations 
are  shown  in  Figure  14. 


FIGURE  14 . MODIFICATIONS  TO  CONTAINER  CONFIGURATION  . 


An  additional  benefit  of  the  fins  is  the  added  static  directional  stability, 
which  should  greatly  aid  the  sirgle-point  suspension  configuration.  A 
test  of  this  configuration  during  the  dynamic  wind  tunnel  tests  showed 
that  the  container  would  remain  in  the  small-end-forward  position  even  when 
started  at  a 90-degree  sideslip.  The  limits  of  stability  for  this 
configuration,  however,  were  not  explored  in  the  tunnel.  It  is  likely 
that  the  addition  of  fins  to  a single-point  suspension  will  increase  the 
limiting  carrying  speed  to  perhaps  60  knots. 
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ACTIVE  STABILITY  AUGMENTATION  USING  CABLE  ACTUATORS 


One  form  of  an  active  stabilization  device  that  was  found  to  be  highly 
effective  is  the  concept  of  driving  the  cables  in  response  to  measured 
load  motions.  In  its  simplest  form,  a cable  guide  is  moved  in  the 
horizontal  plane  by  an  actuator  that  is  driven  by  load  yaw  rate  with  respect 
to  the  helicopter.  This  has  the  effect  of  damping  the  pendulum  and  bifilar 
modes.  If  gains  are  set  to  provide  a damping  ratio  of  0.5  in  hover,  the 
configuration  will  be  easily  handled  in  forward  flight  to  speeds  approaching 
200  knots.  With  a two-point  suspension,  such  a device  would  incorporate 
two  cable  guides  driving  each  cable  laterally  just  below  its  attachment 
point  to  the  helicopter.  The  guides  would  be  driven  laterally  and  collec 
tively  in  response  to  lateral  sway  rate  and  differentially  in  response  to 
load  yaw  rate  (relative  to  helicopter  body  axes).  If  the  gains  are  set 
to  produce  a damping  ratio  of  0.5,  the  limit  cycle  amplitude  is  suppressed 
to  less  than  1/10  of  its  unaugmented  value,  thus  eliminating  any  load 
carrying  speed  limitations.  This  device  was  incorporated  in  the  Model  347 
simulation, and  the  evaluations  revealed  that  speeds  up  to  180  knots  could 
be  reached.  The  maximum  lateral  cable  guide  displacements  called  for  were 
slightly  greater  than  ±2  feet  for  che  transient,  following  turn-on  with  large 
load  motions,  and  less  than  ± foot  for  normal  maneuvering  when  the  system 
was  already  turned  on  prior  to  acceleration  to  high  speed. 

The  advantage  of  this  system  is  its  ability  to  suppress  the  limit  cycle 
for  all  conditions  of  load  angle  of  attack  and  sideslip,  eliminating 
sensitivity  to  maneuvering.  A disadvantage  is,  of  course,  the  cost  and 
complexity  of  the  mechanisms  required  to  sense  the  load  motions  and  actuate 
the  cable  guides. 

MODIFYING  THE  HELICOPTER  STABILITY  AUGMENTATION  SYSTEM 


Modifying  the  helicopter's  stability  augmentation  system  (SAS)  is  an 
attractive  way  of  stabilizing  external  loads.  The  Model  347  helicopter 
system  was  modified  to  include  lateral  SAS  inputs  from  lateral  cable  sway 
rate  and  directional  SAS  inputs  from  the  load  relative  yaw  rate.  No 
shaping  of  these  signals  was  incorporated»  and  the  347 's  normal  low-authority 
SAS  limits  were  retained.  The  gains  were  optimized,  based  on  improving 
the  load  yaw  and  sway  damping  ratios.  Improvements  were  possible  up  to  a 
point, and  then  high-frequency  instabilities  associated  with  reaching  SAS 
limits  were  observed.  The  gains  were  reduced  until  a nominal  damping  ratio 
of  0.3  - 0.4  was  obtained.  This  configuration  was  then  evaluated  and  found 
to  offer  only  slight  increases  in  allowable  carrying  speed. 

A closer  examination  of  the  data  and  pilot  comments  revealed  that  the 
helicopter's  ride  qualities  were  being  compromised  by  the  load  motion  inputs 
into  the  SAS.  Specifically,  the  load  yaw  rates  were  commanding  the 
helicopter  to  "follow"  the  load  in  an  attempt  to  reduce  the  differential 
yaw  motion.  Since  the  pilot's  station  is  about  23  feet  forward  of  the 
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helicopter  center  of  gravity,  this  SAS  activity  showed  up  as  a siie  force 
oscillation  at  the  load  yaw  (bifilar)  frequency.  These  side  forces  became 
uncomfortable  with  increasing  speed  and  caused  the  pilot  to  downrate  nis 
judgment  of  allowable  carrying  speed.  The  lateral  SAS  input  from  load 
sway  rate  did  not  improve  the  situation.  Other  combinations  of  feedbacks 
incorporating  combinations  of  load  position  and  rate  were  attempted, 
without  noticeable  improvements. 

In  a sense,  attempting  to  use  the  relatively  low  bandwidth  helicopter  to 
suppress  the  relatively  high  bandwidth  load  yaw  limit  cycle  motions  does 
not  appear  to  offer  significant  improvements  in  allowable  carrying  speed. 
This  can  be  related  to  the  cable  guide  concept  in  which  the  inertial  cable 
guide  movements  must  be  produced  by  moving  the  whole  helicopter.  These 
movements  were  found  to  be  about  ±1  foot,  which  translates  to  ±5  degrees 
of  helicopter  yaw  angle  for  a cable  separation  of  24  feet. 

ACTIVE  STABILITY  AUGMENTATION  USING  AUXILIARY  AIRFOILS 


The  use  of  an  active  stability  augmentation  system  employing  aerodynamic 
surfaces  driven  by  rate  gyros  is  another  interesting  but  less  practical 
way  of  extending  allowable  carrying  speed.  This  kind  of  device  generally 
employs  a servo-driven  airfoil  surface  placed  behind  the  load  to  introduce 
yawing  moments  proportional  to  yaw  rate.  Because  the  airfoil  must  operate 
in  the  wake  of  the  load,  predicting  its  moment  characteristics  will  be 
difficult.  If  a pure  damping  moment  could  be  produced,  the  device  has  the 
effect  of  augmenting  CN  , which  directly  reduces  the  limit  cycle  amplitude. 

r 

The  amplitude  is  proportional  to  the  reciprocal  of  the  square  root  of  CN  . 

r 

For  the  device  to  be  effective,  must  be  at  least  doubled  from  its 

r 

present  value  of  -0.94.  This  is  a decrement  of  about  -1.  A simple  surface 
that  is  positioned  a container-length  aft  of  the  center  of  mass  may  work; 
however,  this  is  considered  doubtful  because  of  the  highly  separated  and 
turbulent  wake  in  which  the  surface  must  operate.  This  judgment  is  based 
on  drag  chute  test  results  with  the  0.1-scale  container  model  in  the  wind 
tunnel.  If  a drag  device  could  not  stabilize  the  load,  a lifting  device 
such  as  an  airfoil  will  not  do  much  better. 

mmi 

It  is  clear  that  optimizing  the  suspension  system  to  increase  the  bifilar 
frequency,  optimizing  the  angle  of  attack,  and  using  maneuvering  techniques 
can  effectively  extend  carrying  speeds.  Similarly,  it  is  clear  that  the 
addition  of  fins  to  suppress  the  unsteady  aerodynamics  and  the  use  of  more 
exotic  devices  such  as  active  cable  actuators  will  extend  carrying  speeds. 
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Direct  feedback  to  the  onboard  SAS  without  an  extensive  optimization 
effort  produced  only  minor  improvements. 

The  allowable  carrying  speeds  as  determined  from  the  Model  347  simulation 
of  the  empty  8-by-8-20-foot  container  on  a two-point  suspension  are 
listed  in  Table  5.  These  data  are  the  result  of  the  evaluation  of  one 
highly  skilled  pilot  experienced  in  sling  load  flight;  a larger  sample  of 
pilot  opinion  would  enhance  the  results. 
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CONCLUSIONS 


--  * —an  ~ r 


Completion  of  the  basic  study  invo’ving  flight  simulation  of  the  Model 
347  tandem-rotor  helicopter  as  well  ,'s  the  closely  related  simulation 
studies  generated  data  which  support  the  following  conclusions: 


1.  The  Model  347  flight  simulation  program  produced  an  effective 
representation  of  a large  tandem-rotor  helicopter. 

2.  Unsteady  aerodynamic  effects  are  a predominant  factor  in  helicopter 
sling  load  dynamics. 

3.  Use  of  a nominal  ICdegree  load  angle  of  attack  in  pitch  and  yaw 
will  minimize  the  unsteady  effects  of  sling  loads. 

4.  Sling  load  dynamics  can  also  be  minimized  by  using  wide  cable 
separation  and  short  cable  lengths. 
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APPENDIX 


ADVANCED  TANDEM-ROTOR  HELICOPTER  HYBRID  COMPUTER  MODEL 


The  schematic  diagrams  for  the  full  flight  envelope  mathematical  model 
of  the  basic  347  aircraft  are  presented  in  this  appendix.  Functional 
block  diagrams  included  at  the  beginning  of  this  section  present  a genera- 
lized picture  of  how  the  various  portions  of  the  mathematical  model  are 
related.  The  relationship  of  the  various  digital  and  analog  computers 
and  where  the  various  mathematical  computations  were  assigned  are  also 
shown.  The  longitudinal,  lateral,  and  directional  Stability  Augmentation 
System,  rotor  speed  governor,  and  engine  dynamics  and  stick  force  models 
are  also  shown  in  functional  block  diagram  format. 

Detailed  diagrams  showing  appropriate  analog  computer  component  patching 
follow.  The  equations  which  are  being  modelled  have  been  annotated  on 
the  appropriate  drawings. 

A complete  list  of  symbols  for  the  report  body  and  for  the  appendix  is 
included  following  the  diagrams. 
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AT  (OKI 


Symbol 


Axis  transformation  matrix  from  inertial  to 
aircraft  body  axes 

Transpose  of  the  orthogonal  matrix  of  direction 
cosine  transformation  from  inertial  to  aircraft 
body  axes 

Lateral  cyclic  pitch  in  body  axis  * 

Axis  transformation  from  inertial  to  load 
body  axes 

Transpose  of  the  orthogonal  matrix  of  direction 
cosine  transformation  from  inertial  to  load  tody 


DYNAMIC 


Longitudinal  cyclic  pitch  in  body  axis 
Sling  load  unsteady  yawing  moment  coefficient 


DAMPING 


Sling  load  yawing  moment  coefficient 

Total  aerodynamic  sling  load  yawing  moment 
coefficient  due  to  static,  ■ .mic  and  yaw 
rate  components 

Sling  load  yawing  moment  coefficient  due 
to  sideslip 

Sling  load  yawing  moment  coefficient  due  to 

yaw  rate’  C ib. 


Sling  load  yawing  moment  coefficient  due  to  non- 
dimensional  yaw  rate,  rb^ 


STATIC 


Sling  load  static  load  yawing  moment  coefficient 


Symbol 

Description 

1MI 

°Y 

Sling  load  side  force  coefficient  due  to 

1/Rad 

yb 

sideslip 

D/q 

Nondimensional  fuselage  drag  force  in 
body  axis 

Ft2 

°LT 

Sling  load  drag  in  body  axis  of  load 
<Dlt  » D/q) 

Ft2 

G(  ju>) 

Amplitude  ratio 

Db 

V «R 

Rotor  (H)  drag  force  in  S.N.P.  wind  axis 

Lb 

Hz 

Frequency  in  cycles  per  second 

- 

I 

XX 

Airframe  moment  of  inertia  about  X axis 

Slug-Ft 

I 

yy 

Airfiame  moment  of  inertia  about  Y axis 

Slug-Ft 

i 

zz 

Airframe  moment  of  inertia  about  Z axis 

Slug-Ft 

I 

xz 

Airframe  product  of  inertia  about  XZ  plane 

Slug-Ft 

lZ  °r  JL 

Yaw  moment  of  inertia  of  sling  load 

Slug-Ft 

n n 

Linear  nrt  cable  stretch  (and  angular  rate) 
constant 

Lb/Ft 

Lb/Ft/Sec 

L 

Total  airframe  rolling  moment  in  aircraft 
body  axis 

Ft-Lb 

lat 

n 

nch  cable  tension  rolling  moment  in  helicopter 
axes 

Ft-Lb 

lbt 

n 

nth  cable  tension  rolling  moment  in  load 
body  axes 

Ft-Lb 

M 

Total  airframe  pitching  moment  in  aircraft 
body  axis 

Ft-Lb 

Mirr, 

AT 

n 

nth  cable  tension  pitching  moment  in  helicopter 
axes 

Ft-Lb 
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Symbol 

Description 

Unit 

| ^BT 

n 

nt!l  cable  tension  pitching  moment  in  load 
body  axes 

Ft-Lb 

N 

Total  airframe  yawing  moment  in  aircraft 
body  axis 

Ft-Lb 

N/q 

Non dimensional  fuselage  yawing  moment 
in  body  axis 

Ft3 

naero 

Total  sling  load  aerodynamic  yawing  moment 
in  body  axes  of  the  load 

Ft-Lb 

1 nat 

n 

nth  cable  tension  yawing  moment  in  helicopter 
axes 

Ft-Lb 

| nbt 

n 

nth  cable  tension  yawing  moment  in  load 
body  axes 

Ft-Lb 

ngables 

Sling  load  cable  yawing  moment 

Ft-Lb 

i 

N 

LT 

Sling  load  yawing  moment  in  body  axis  of 
load  (Nlt  = N/q) 

Ft3 

N 

TOTAL 

Total  sling  load  aerodynamic  and  cable 
yawing  moment  in  body  axes  of  the  load 

Ft-Lb 

1 

nunsteady 

Sling  load  unsteady  yawing  moment 

Ft-Lb 

p,p 

Component  of  helicopter  angular  rate  (or 
acceleration)  about  aircraft  X axis 

Rad/Sec 

Rad/Sec 

Q,Q 

Component  of  helicopter  angular  rate  (or 
acceleration)  about  aircraft  Y axis 

Rad/Sec 

Rad/Sec 

qf’\ 

Rotor  torque  required 

Ft-Lb 

• 

R,R 

Component  of  helicopter  angular  rate 
(or  acceleration)  about  aircraft  Z axis 

Rad/ Sec 
Rad/Sec 

SL 

Load  container  planform  area,  x b^ 

Ft2 

ST 

Strouhal  number, 

Rad 

T 

0 

Cable  tension  force 

Lb 
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Symbol 


IML 


tf’tr 


Rotor  thrust  in  S.N.P.  wind  axis 


Tension  force  of  the  nth  cable 


VT  °r  Vc 


WL  or  W 


Total  aerodynamic  velocity 

Total  steady  wind  velocity 

Weight  of  sling  load 

Total  force  along  helicopter  X axis 


Helicopter  body  axis  coordinate  in  X 
direction 

nth  cable  tension  force  along  helicopter 
X axis 

Load  body  axis  coordinate  in  X direction 

nth  cable  tension  force  along  load  x axis 

Load  container  half-length 

Horizontal  distance  between  cable  and  load 
center  of  mass 

Total  force  along  helicopter  Y axis 

Nondimensional  fuselage  side  force  in 
body  axis 

Helicopter  body  axis  coordinate  in  y direction 

nch  cable  tension  force  along  helicopter 
y axis 

Load  body  axis  coordinate  in  y direction 
nth  cable  tension  force  along  load  y axis 


YF’YR 


Rotor  side  force  in  S.N.P.  wind  axis 


Symbol 

PcacrlptlflB 

Unit 

ylt 

Sling  load  side  force  iu  body  axis  of  load 

(VLI  » v/q> 

Ft2 

Z 

Total  force  along  helicopter  Z axis 

Lb 

ZA 

Helicopter  body  axis  coordinate  in  Z 
direction 

- 

ZAT 

n 

nth  cable  tension  force  along  helicopter 
Z axis 

Lb 

ZB 

Load  body  axis  coordinate  in  Z direction 

- 

ZBT 

n 

nth  cable  tension  force  along  load  Z axis 

Lb 

a ,a  .a 
° °F  °R 

Rotor  coning  angle 

Rad 

ai*a»  »ai 

1F  R 

Longitudinal  first  harmonic  flapping 
coefficient 

Rad 

Ybi  »bi 

F R 

Lateral  first  harmonic  flapping  coefficient 

Rad 

bL 

Sling  load  container  width 

Ft 

dt 

Rate  of  change  of  time 

Sec 

f 

Sling  load  unsteady  yawing  moment  function 

- 

F 

Acceleration  due  to  gravity 

Ft/Sec 

j 

V^r 

- 

l,ln 

Straight  line  distance  between  upper  cable 
(or  nth  cable)  and  lower  cable  (or  nch  cable) 
sling  attach  points 

Ft 

Bifilar  length 

Ft 

*l 

Sling  load  container  length 

Ft 

*n 

Rate  of  change  of  straight  line  distance 
between  upper  nth  cable  and  lower  nth  cable 
sling  attach  points 

Ft/Sec 
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Symbol 

Description 

Unit 

l ,1 

o’  o 

n 

Zero-tension  straight  line  distance 
between  upper  cable  (or  rf*  cable)  and 
lower  cable  (or  nch  cable)  sling  attach  points 

Ft 

1 P 

Pendulum  length 

Ft 

m 

Mass  of  helicopter 

Slugs 

PA 

Component  of  helicopter  angular  rate  about 
helicopter  X axis 

Rad/Sec 

PB 

Component  of  load  angular  rate  about  load 
X axis 

Rad/Sec 

“A 

Component  of  helicopter  angular  rate  about 
helicopter  y axis 

Rad/Sec 

Component  of  load  angular  rate  about  load 
y axis 

Rad/Sec 

rA 

Component  of  helicopter  angular  rate  about 
helicopter  Z axis 

Rad/Sec 

rB 

Component  of  load  angular  rate  about  load 
Z axis 

Rad/Sec 

t 

Time 

Sec 

• 

u,  u 

u or  u. 
A 

Linear  inertial  velocity  (and  acceleration) 
components  along  helicopter  X axis 

Ft/Sec, 

Ft/Sec 

UB 

Linear  inertial  velocity  along  load  X axis 

Ft/Sec 

• 

u„.  ,xni 
BA  * BA 
n n 

Difference  between  (nth  cable)  lower  attach 
point  velocity  and  (nth  cable)  upper  attach 
point  velocity,  along  helicopter  X axis, 
relative  to  helicopter  axes 

Ft/Sec 

• 

V,  V 

v or  v. 

A 

Linear  inertial  velocity  (ai.  acceleration) 
components  along  helicopter  Y axis 

Ft/Sec„ 

Ft/Sec 

VB 

Linear  inertial  velocity  along  load  y txis 

Ft/Sec 

• 

VBA  ,yBA 
n n 

Difference  between  (nth  cable)  lower  attach 
point  velocity  and  (nth  cable)  upper  attach 
point  velocity,  along  helicopter  y axis, 
relative  to  helicopter  axes 

Ft/Sec 
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Symbol 

• 

w,  w 
w or  vA 

WB 

WBA  »ZBA 
n n 


x 


XA*XA 


y 


yA*yA 


yB*yB 


z 


Description 

Linear  inertial  velocity  (and  acceleration) 
components  along  helicopter  Z axis 

Linear  inertial  velocity  along  load  Z axis 

Difference  between  (nch  cable)  lower  attach 
point  velocity  and  (n:h  cable)  upper  attach 
point  velocity,  along  helicopter  Z axis, 
relative  to  helicopter  axes 

Body  axis  coordinate  in  X direction 

Cable  (or  nth  cable)  upper  attach  point 
along  helicopter  X axis  relative  to  the 
helicopter-fixed  axes 

Cable  (or  nch  cable)  lower  attach  point  along 
load  X axis  relative  to  the  load-fixed  axes 

Difference  between  nch  cauls  lower  attach 
point  position  and  nth  cable  upper  attach 
point  position,  along  helicopter  x axis, 
relative  to  helicopter  axes 

Body  axis  coordinate  in  y direction 

Cable  (or  nth  cable)  upper  attach  point  along 
helicopter  y axis  relative  to  the  helicopter  - 
fixed  axes 

Cable  (or  nth  cable)  lower  attach  point  along 
load  y axis  relative  to  the  load-fixed  axes 

Difference  between  nth  cable  lower  attach 
point  position  and  nth  cable  upper  attach 
point  position,  along  helicopter  y axis, 
relative  to  helicopter  axes 

Body  axis  coordinate  in  z direction 

Cable  (or  nch  cable)  upper  attach  point  along 
helicopter  z axis  relative  to  the  helicopter  - 
fixed  axes 

Lower  cable  (or  nth  cable)  lower  attach  point 
along  load  z axis  relative  to  the  load-fixed 
axes 


Unit 

Ft/Sec2 

Ft/Sec 

Ft/Sec 

Ft/Sec 


Ft 


Ft 
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Symbol 

Pescrijtioa 

Unit 

ZBA 

n 

Difference  between  ntk  cable  lower  attach 
point  position  and  nch  cable  upper  attach 
point  position,  along  helicopter  z axis, 
relative  to  helicopter  axes 

Ft 

$ G(juj) 

Phase  angle 

Deg 

ttFUS 

Angle  of  attack  of  fuselage  (includes  rotor 
downwash  effects) 

Rad 

aL»®L 

Sling  load  angle  of  attack  (and  angular 
rate) 

Rad 

Rad/Sec 

9 

Fuselage  sideslip  angle 

Rad 

3' 

Local  sideslip  angle  of  rotors 

Rad 

Ml 

p 

'command 

^COMMAND 

Sling  load  angle  (and  angular  rate)  of  slip 
Sling  load  sideslip  (and  angular  rate)  command 

Rad 

Rad/Sec 

Rad 

Rad/Sec 

3 ,P 

€ € 

Sling  load  sideslip  (and  angular  rate)  error 
between  command  and  actual  sideslip  angle 
(and  angular  rate) 

Rad 

Rad/Sec 

63 

Forward  rotor  delta  three  designation 

Deg 

6b 

Longitudinal  stick  cockpit  control 

In 

6c 

Collective  stick  cockpit  control 

In 

6r 

Directional  pedal  cockpit  control 

In 

6S 

Lateral  stick  cockpit  contro1 

In 

e 

Euler  pitch  angle 

Rad 

e ,e 

°F  °R 

Root  collective  pitch  (0  is  corrected 

for  53)  °F 

Rad 

0oFM,eoRM 

"Rotor  map"  757„  blade  radius  collective 
pi  tch 

Rad 

Symbol 

Deserf ption 

Unit 

Complete  inflov.  ratio 

- 

MT 

Transpose  of  the  orthogonal  matrix  of 
direction  cosines  of  the  taut  cable  after 
bein';  rotated  through  the  angles  and 

X2n 

H’Xn 

Longitudinal  cable  (or  nth  cable)  angle 
parallel  to  the  XZ  plane  of  the  helicopter 

Rad 

Longitudinal  nth  cable  velocity  parallel 
to  the  XZ  plane  of  the  helicopter 

Rad/Sec 

X2»Xn 

Lateral  cable  (or  nth  cable)  angle  parallel 
to  the  YZ  plane  of  the  helicopter 

Rad 

X2n 

Lateral  nth  cable  velocity  parallel  to  the  YZ 
plane  of  the  helicopter 

Rad/Sec 

H 

Advance  ratio  u , u 

r ft 

- 

QR  QR 

0 

Euler  roll  angle  and  angular  rate 

Rad 

VV*'l 

Sling  load  yaw  angle  (and  yaw  rate) 
(and  yaw  acceleration) 

Rad 

Rad/Sec 

Rad/Sec 

ylimit 

Allowable  sling  load  limit  cycle  amplitude 

Rad 

^max 

Maximum  amplitude  of  the  sling  load  yaw 
limit  cycle 

Rad 

Y 

w 

Direction  from  which  steady  wind 
blows  expressed  from  0 to  ± 180° 

Rad 

(JU 

Frequency 

Rad/Sec 

("Y 

Bifilar  frequency 

Rad,  Sec 

ct 

Real  part  of  complex  root 

- 
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Symbol 

Description 

iiaJLL 

Aic 

Lateral  cyclic  pitch  in  body  axis 

Rad 

AICF' 

Lateral  cyclic  pitch  control  input  in 
body  axis  - forward  rotor 

Rad 

A-  _ _ 

ICR 

Lateral  cyclic  pitch  control  input  in 
body  axis  - aft  rotor 

Rad 

aicfr» 

Lateral  cyclic  pitch  in  the  shaft  normal 
plane  wind  axis  - forward  rotor 

Rad 

aicfr 

Lateral  cyclic  pitch  in  the  shaft  normal 
plane  wind  axis  corrected  for  63  hinging 
and  control  phasing  - forward  rotor 

Rad 

aicrr 

Lateral  cyclic  pitch  in  the 
plane  wind  axis  - aft  rotor 

Rad 

asl 

Sling  load  coefficient 

- 

bic 

Longitudinal  cyclic  pitch  in  body  axis 

Rad 

bicf» 

Longitudinal  cyclic  pitch  control  input 
in  body  axis  - forward  rotor 

Rad 

bicr 

Longitudinal  cyclic  pitch  control  input 
in  body  axis  - aft  rotor 

Rad 

BICFR' 

Longitudinal  cyclic  pitch  in  the  shaft 
normal  plane  wind  axis  - aft  rotor 

Rad 

bicfr 

Longitudinal  cyclic  pitch  in  the  shaft 
normal  plane  wind  axis  corrected  for  4^ 
hinging  and  control  phasing  - aft  rotor 

Rad 

bicrr 

Longitudinal  cyclic  pitch  in  the  shaft 
normal  plane  wind  axis  - aft  rotor 

Rad 

bicfrm 

"Rotor  map"  shaft  normal  plane  wind  axis 
longitudinal  cyclic  pitch  corrected  for 
6 ^ hinging  and  control  phasing  - forward  rotor 

Rad 

bicrrm 

"Rotor  map"  shaft  normal  plane  wind  axis 
longitudinal  cyclic  px^-ch  - aft  rotor 

(bicrrm  = bicrr) 

Rad 
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feata  JfcfiiiMlMi 


Svmbol 

Description 

Unit 

bitf 

Longitudinal  cyclic  trim  in  the  hody 
axis  - forward  rotor 

Deg 

bitr 

Longitudinal  cyclic  trim  in  the  body  axis  ~ 
aft  rotor 

Deg 

bsl 

Sling  load  coefficient 

- 

chf,chr 

Rotor  H (drag)  force  coefficient  in  shaft 
normal  plane  wind  axis 

- 

CH  /a 

"Rotor  map"  derived  rotor  drag  force  coefficient 
in  shaft  normal  plane  wind  axes 

- 

CHF/a 

"Rotor  map"  derived  rotor  drag  force  coefficient 
in  shaft  normal  plane  wind  axis 

- 

HR/a 

2CHF/aa 

2C 

Z HR/ac 

Rotor  H force  coefficient  calculated  with 
classical  equation  - shaft  normal  plane 
wind  axis 

— 

hftot 

"Rotor  map"  derived  rotor  drag  force  coefficient 
corrected  for  shaft  pitching  angular  velocity 

- 

hrtot 

CPF’CPR 

"Rotor  map"  power  coefficient,  adjusted  for 
compressibility(Mach  number  correction  applied) 

- 

CPF/a 

r> 

Nondimensional  "rotor  map"  power  coefficient 
adjusted  for  Mach  number  correction 

- 

PR/a 

CPFl/a 

r 

"Rotor  map"  nondimensional  power  coefficient 
(uncorrected) 

- 

PRl/a 

CQ  ,CQ 
yF  yR 

Rotor  torque  coefficient  (C^  = Cp) 

• 

Rotor  aero  torque  required  coefficient 
| computed  with  "classical"  equation  - forward 

F rotor 
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Symbol 


Pescrigtifln 


Unit. 


Rotor  aero  torque  required  coefficient 
computed  with  "classical”  equation  - 
aft  rotor 


T 

R/c 

CTF1/ct 

CTRl/a 


aCT 


a a 


Sling  load  coefficient 

Rotor  thrust  coefficient  in  shaft  normal 
plane  wind  axis 

"Rotor  map"  thrust  coefficient  in  shaft  normal 
plane  wind  axis 

"Rotor  map"  thrust  coefficient  resolved 
in  .o  shaft  normal  plane  wind  axis 


"Rotor  map"  thrust  coefficient  as  stored 
in  map  tables  (rotor  wind  axis) 


Rotor  thrust  coefficient  computed  with 
classical  equation  (shaft  normal  plane 
wind  axis) 


Rotor  propulsive  force  coefficient  in 
rotor  wind  axis 


"Rotor  map"  propulsive  force  coefficient  as 
stored  in  map  tables  (rotor  wind  axis) 


Rotor  side  force  coefficient  in  shaft 
normal  plane  wind  axis 


2Cy  Rotor  side  force  coefficient  computed 

F with  "classical"  equation  (shaft  normal 

aa  plane  wind  axis) 


ac 
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Symbol 

D/q 


FFR 


FRF 


LT 


SL 


SL 


F 

F. 


JBR 


JFR 


SL 

gi’gq 

gj8 


AP 


Description 

Nondimen sional  fuselage  drag  force  in 
body  axis 

Forward  on  aft  rotor  interference  factor 
(percent  of  fwd  rotor  uniform  induced 
velocity  seen  by  aft  rotor) 

Aft  on  forward  rotor  interference  factor 
(percent  of  aft  rotor  induced  velocity 
seen  by  fwd  rotor ) 

Sling  load  drag  in  body  axis  of  fuselage 

Sling  load  drag  in  body  axis  of  load 

(Dlt  = D/q) 

Sling  load  coefficient 

Sling  load  coefficient 

4 

Rotor  force  coefficient  parameter  (prrRg  ) 
Longitudinal  stick  breakout  force 
Longitudinal  stick  frictional  force 


Halt 


Ft 


Lb/Sec2 

Lb 

Lb 


Sling  load  coefficient 

Engine  governor  torque/rpm  gain  Ft-Lb/Rad/Sec 


Sideslip  SAS  gain  - 347  analog  control  In. Eauiv. Pedal 

system  Rad  Sideslip 


Sideslip  SAS  gain  as  mechanized  in 
347  Demonstrator  Task  III  control 


In. Eauiv. Pedal 
Lb/Ft2  of  AP 


Sling  load  coefficient 


SAS  gain 
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Symbol 

H 


hf,Hr 


SL 


V1. 


I 


SL 


I 


T 


I 

xx 


I 

yy 


i 

xz 


J 


L 


K 


K 

kdf,kdr 


£,  L 


Bfififcatttifln. 

Pressure  altitude 

Rotor  (H)  drag  force  in  shaft  normal 
plane  wind  axis 

Sling  load  coefficient 

Individual  rotor  and  drive  system 
polar  moment  of  inertia 

Sling  load  coefficient 

Engine  power  turbine  moment  of  inertia 

Airframe  moment  of  inertia  about  X axis 

Airframe  moment  of  inertia  about  Y axis 

Airframe  moment  of  inertia  about  Z axis 

Airframe  product  of  inertia  about  X plane 

Yaw  moment  of  inertia  of  sling  load 

Engine  power  turbine  internal  damping 
constant 

Factor  in  sling  load  equations  of  motion 

Special  error  rate  damping  feedback  constant 
for  rotor  speed  governor  stabilization 

Rotor  shafting  spring  constant  in  rotor 
speed  governor  circuit 

Sling  load  coefficient 

(-)  Tangent  of  forward  rotor  delta  three 
ingle 

Total  airframe  rolling  moment  in  aircraft 
body  axis 


Unit 

Ft 

Lb 


Slug-Ft2 

Slug-Ft2 

Slug-Ft2 

Slug-Ft2 

Slug-Ft2 

Slug-Ft2 

Slug-Ft2 

Ft-Lb/Rad/Sec 

Ft-Lb/Rad/Sec 

Ft-Lb 


Ft-Lb 
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Symbol 

Ppsrri-Dtion 

Nondimens ional  fuselage  rolling  moment, 
body  axis 

L/q 

Nondimensional  fuselage  lift  in  body  axis 

L' 

Rolling  moment  contribution  of  rotor 
T,rl,  and  Y forces  about  aircraft  c.g. 

lfus 

Fuselage  rolling  moment  in  body  axis 

lhf’lhr 

Rotor  lateral  hub  moment  due  to  flapping 
expressed  in  shaft  normal  nlane  wind  axis 

ll 

Average  sling  load  cable  length  below 
attachment  point 

M 

Total  airframe  pitching  moment  in  aircraft 
body  axis 

M/q 

Nondimensional  fuselage  pitching  moment 
in  body  axis 

M' 

Pitching  moment  contribution  of  rotor  T, 
H,  and  Y forces  about  aircraft  c.g. 

M 

FUS 

mh 

Fuselage  pitching  moment  in  body  axis 

Helicopter  mass  in  sling  load  equations 
of  motion 

^F^HR 

Rotor  longitudinal  hub  moment  due  to 
flapping  - expressed  in  shaft  norma 
plane  wind  axis 

MoF’MoR 

Advancing  tip  Mach  number  of  rotor  for 
"rotor  map"  calculations 

^REFp’ 

Reference  advancing  tip  Mach  number  of 
rotor  at  which  maps  were  calculated 

MREFr 

msl 

Sling  Load  coefficient 

UnU 


Ft' 


Ft 

Ft-Lb 


Ft-Lb 

Ft-Lb 


Ft 

Ft-Lb 

Ft3 

Ft-Lb 

Ft-Lb 

Slugs 

Ft-Lb 


M 

M 
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Symbol 

Description 

Unit 

*r 

Sling  load  coefficient 

- 

M 

w 

Mass  moment  of  blade  about  flapping 
hinge  used  in  hub  moment  calculations 

Slug-Ft 

N 

Total  airframe  yawing  moment  in  aircraft 
body  axis 

Ft-Lb 

N/q 

Nondimens ional  fuselage  yawing  moment 
in  body  axis 

Ft3 

N' 

Yawing  moment  contribution  of  rotor  T, 
H,  and  Y forces  about  aircraft  c.g. 

Ft-Lb 

nfus 

Fuselage  yawing  moment  in  body  axis 

Ft-Lb 

nl 

Sling  load  yawing  moment  resolved 
into  fuselage  body  axis 

Ft-Lb 

nlt 

Sling  load  yawing  moment  in  body  axis 
of  load  (Nlt  = N/q) 

-j 

Ft” 

nsl  < 

Sling  load  coefficient 

- 

P’PA’P 

Component  of  helicopter  angular  rate 
(or  acceleration)  about  aircraft  X 

Rad/Sec 

axis 

Rad/Sec 

PSL 

Sling  load  coefficient 

- 

Q,QA,p 

Component  of  helicopter  angular  rate 
(or  acceleration)  about  aircraft  Y 

Rad/Sec 

axis 

Rad/Sec 

^AERO  ’ 
r 

Rotor  aerodynamic  torque  required 

Ft-Lb 

QAER0o 

K 

% 

Total  initial  gas  generator  output  torque 

Ft-Lb 

Q’g 

Engine  output  torque 

Ft-Lb 

QG  ’ V 

Rotor  shafting  windup  torque 

Ft-Lb 
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Symbol 

QSL 

r,  rar 


S 

SSL 

tf’tr 


v (nRj) 


v 

w 


V 

x 


V 

y 


wT 


WSL 
X,  X 


Description 

Sling  load  coefficient 

Component  of  helicopter  angular 
rate  (or  acceleration)  about  aircraft 
Z axis 

Rotor  blade  radius 

Vertical  distance  between  sling  load 
hook  attachment  and  aircraft  c.g. 

Laplace  operator 

Sling  load  coefficient 

Rotor  thrust  in  shaft  normal  plane 
wind  axis 

Absolute  temperature 
Sling  load  coefficient 


Average  rotor  induced  velocity  used 
in  computing  rotor  downwash  on 
fuselage 

Total  aerodynamic  velocity 
Total  steady  wind  velocity  t 

Longitudinal  aircraft  velocity  in  earth 
coordinates  resolved  along  the  helicopter 
x axis 

Lateral  aircraft  velocity  in  earth  coordinates 
resolved  along  the  helicopter  y axis 

Weight  of  sling  load 

Sling  load  coefficient 

Total  force  (and  linear  velocity)  Lb 

along  helicopter  X axis 


Unit 

Rad/Sec_ 

Rad/Sec 

It 

Ft 

Lb 

°F 

Ft/Sec 

Ft/Sec 

Ft/Sec 

Ft/Sec 

Ft/Sec 

Lb 

Ft/Sec 
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Symbol 

XFUS 

^2 

XSL 

Y.i 

Vyr 


Pessriptign 

Fuselage  X force  expressed  in  body  axis 

Simulator  fore  and  aft  position 

Sling  load  coefficient 

Total  force  (and  linear  velocity) 
along  helicopter  Y axis 

Rotor  side  force  in  shaft  normal  plane 
wind  axis 


Y 

FUS 


Fuselage  Y force  expressed  in  body  axis 


Y/q 


Nondimens ional  fuselage  side  force 
in  body  axis 


Sling  load  side  force  in  body  axis  of 
fuselage 


Sling  load  side  force  in  body  axis 
of  load  (Ylt  = Y/q) 

Simulator  side  position 

Total  force  (and  linear  velocity) 
along  helicopter  Z axis 


ZFUS 


M2 


a 


Fuselage  Z force  expressed  in  body  axis 

Simulator  vertical  position 

Rotor  blade  section  lift  curve  slope 

Cable  separation  (distance  from  forward 
to  aft  cable  attachment) 


a ,a  Rotor  coning  angle 

°F  °R 

a.  ,a.  Longitudinal  first  harmonic  flapping 

F R coefficient 


b 


Number  of  blades  per  rotor 


Unit 

Lb 

Ft 


Lb,  Ft/Sec 
Lb 
Lb 


Lb 


Ft 

Lb,  Ft/Sec 


Lb 

Ft 

Rad 

Ft 


Rad 

Rad 


blF,blR 


Lateral  first  harmonic  flapping  coefficient 


Rad 


Symbol 


Description 


Halt 


c 

e 

f 

i 

g 

h 


Vr 


h 

c 


1F*lR 

V^R 


q 


u,  u 


UF’UR 


Blade  chord 

Rotor  flapping  hinge  offset 

Equivalent  flat  plate  area  (drag) 
of  fuselage 

Acceleration  due  to  gravity 

Vertical  rate  of  climb  in  inertial 
axis  (h  = -Z) 

Vertical  distance  from  helicopter  c.g. 
to  rotor  hub  measured  parallel  to 
helicopter  Z axis 

Vertical  distance  between  the  required 
helicopter  c.g.  and  the  wind  tunnel 
model  test  c.g. 

Rotor  shaft  incidence  angles 

Longitudinal  distance  from  helicopter 
c.g.  to  projection  of  rotor  hub  on  the 
X axis 

Longitudinal  distance  between  the  required 
helicopter  c.g.  and  the  wind  tunnel  test 

c.g. 


Mass  of  sling  load 

Local  dynamic  pressure  on  fuselage 

Local  dynamic  pressure  as  seen  by  sling 
load 

Linear  inertial  velocity  (and 
acceleration)  components  along  helicopter 
X axis 

Longitudinal  velocity  component  at 
the  rotor  parallel  to  local  wind  and 
in  shaft  normal  plane 


Ft 

Ft 

Ft2 

Ft/Sec2 

Ft/Sec 

Ft 

Ft 


Deg 

Ft 

Ft 


Slugs 

Lb/Ft2 

Lb/Ft2 

F c/Sec 

Ft/Sec2 

Ft/Sec 
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Svmbol 


V\ 


UF  * °R 
F2  R2 


u ,u 
gF  gR 


Description 

Longitudinal  velocity  component  in 
body  axis  corrected  for  rotor  angular 
pitching  rate  (includes  wind) 

Longitudinal  velocity  component  at 
rotor  corrected  for  rotor  angular  rate 
and  resolved  into  shaft  normal  plane 

’•1st  velocity  component  parallel  to 
!i_  •'•raft  X axis  as  seen  by  forward 
rotor  or  aft  rotor 


u^  Sling  load  longitudinal  velocity  expressed 

in  fuselage  body  axis  system 

UpQ  Longitudinal  specific  force 

iipop  Filtered  specific  force 

u Simulator  longitudinal  specific  force 

F°FB  feedback 


u 


T 


u 

V 


v,v 


v 


F 


2 


Total  fuselage  aerodynamic  velocity 
component  along  body  X axis  (includes 
wind) 

Component  of  steady  wind  velocity 
along  aircraft  X axis 

Linear  inertial  velocity  (and  acceleration) 
components  along  helicopter  Y axis 

Lateral  velocity  in  body  axis 
corrected  for  rotor  angular  rolling 
and  yawing  rates  (includes  wind) 

Lateral  velocity  component  at  rotor 
corrected  for  rotor  angular  rate 
and  resolved  into  shaft  normal  plane 

Gust  velocity  component  parallel  to 
aircraft  Y axis  as  seen  by  forward  or  aft 
rotor 


Mil 

Ft/Sec 

Ft/Sec 

Ft/Sec 

Ft/Sec 

Ft/Sec2 

Ft/Sec2 

Ft/Sec2 

Ft/Sec 

Ft/Sec 

Ft/Sec~ 

Ft/Sec 

Ft/Sec 

Ft/Sec 

Ft/Sec 
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M 





Svmbol 

Description 

JMt 

V1NDUCED 

Rotor  uniform  induced  velocity 

v J?. i nn 

v2  n2+x2 

Ft/Sec 

VL 

Sling  load  lateral  velocity  expressed 
in  fuselage  body  axis  system 

Ft/Sec 

• 

V 

po 

Lateral  specific  force 

Ft/Sec 

VpoF 

Filtered  specific  force 

Ft/Sec 

• 

V 

P°FB 

Simulator  lateral  specific  force  feedback 

Ft/Sec 

VT 

Total  fuselage  aerodynamic  velocity 
component  along  body  Y axis  (includes 
wind) 

Ft/Sec 

V 

w 

Component  of  steady  wind  velocity 
along  aircraft  Y axis 

Ft/Sec 

w,w 

Linear  inertial  velocity  (and  acceleration) 
components  along  helicopter  Z axis 

Ft/Sec 

VWR 

WF  = WF2’  WR  = WR2 

Ft/Sec 

WF  WR 

Fl,  R1 

Vertical  velocity  in  body  axis  corrected 
for  rotor  angular  pitching  rate 

Ft/Sec 

WF  ’ WR 
2 2 

Vertical  velocity  component  at  rotor 
corrected  for  rotor  angular  rate  and 
resolved  into  shaft  normal  plane 

Ft/Sec 

w ,w 
gF  gR 

Gust  velocity  component  parallel  to 
aircraft  Z axis  as  seen  by  forward 
or  aft  rotor 

Ft/Sec 

Sling  load  vertical  velocity  expressed 
in  fuselage  body  axis  system 

Ft/Sec 

• 

W_ 

Po 

Normal  specific  force 

Ft/Sec 

"Po. 

Filtered  normal  specific  force 

Ft/Sec 
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Symbol 

Description 

Unit 

• 

P°FB 

WT 

Simulator  normal  specific  force  feedback 

Total  fuselage  aerodynamic  velocity 
component  along  body  Z axis  (includes 
wind) 

Ft/Sec 

Ft/Sec 

WT* 

Total  vertical  velocity  acting  on  the 
fuselage  modified  for  rotor  downwash 

Ft/Sec 

w 

w 

Component  of  steady  wind  velocity 
along  aircraft  Z axis 

Ft/Sec 

X 

Body  axis  coordinate  in  X direction 

- 

XF,XR 

Rotor  shaft  normal  plane  wind  axis 
coordinate  in  X direction 

y 

Body  axis  coordinate  in  Y direction 

- 

yF’yR 

Rotor  shaft  normal  plane  wind  axis 
coordinate  in  Y direction 

- 

2 

Body  axis  coordinate  in  Z direction 

- 

ZF’ZR 

Rotor  shaft  normal  plane  wind  axis 
coordinate  in  Z direction 

- 

aCAF’aCAR 

Control  axis  angle  of  attack  in 
"rotor  map"  calculations 

Rad 

aFUS 

Angle  of  attack  of  fuselage  (includes 
rotor  downwash  effects) 

Rad 

al 

Sling  load  angle  of  attack 

Rad 

aSF’aSR 

Angle  of  attack  of  the  rotor  shaft 
normal  plane 

Rad 

^FUS 

Fuselage  sideslip  angle 

Rad 

Local  sideslip  angle  of  forward  and 
aft  rotors 

Rad 
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Unix 


Symbol 

V 


6b,6b 

p 

6b 

DASH 


Pefifiiijtlpa 

Sling  load  angle  of  slip 

Rotor  blade  Lock  number  (blade  inertia 
parameter) 

Total  left  or  right  lateral  swing 
angle  of  the  forward  sling  load 
cable  as  viewed  from  rear  of  aircraft 

Total  left  or  right  lateral  swing 
angle  of  the  aft  sling  load  cable 
as  viewed  from  rear  of  aircraft 

Average  rotor  blade  section  prorile 

? 

draft  coefficient  (6  = 6q  + 6^  a ) 

Constant  term  in  section  drag  expression 

Coefficient  of  pressure  drag  due  to 
lift  term  in  section  drag  expression 

Forward  rotor  delta  three  designation 

Longitudinal  stick  (and  pilot)  control 
input 

Longitudinal  dash  control  input 


Longitudinal  stick  driver  command  input 


Longitudinal  SAS  control  input 


Longitudinal  stick  trim  control  input 


Collective  stick  (and  pilot)  control 
input 


Longitudinal  stick  magnetic  brake  control 
input 


Rad 


Rad 


Rad 


Deg 


In. 


In. 


In. 


In.  of 
equivalent 
cockpit 
control 


In . 


In . 


In . 
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Symbol 


Description 


Rudder  pedal  magnetic  brake  control  input 


6r»6r 

P 


Lateral  stick  magnetic  brake  control  input 

Directional  pedal  (and  pilot)  control 
input 

Rudder  pedal  driver  command  input 
Directional  SAS  control  input 


6 


*T 


6S’6S 

P 


Rudder  pedal  trim  control  input 
Lateral  stick  (and  pilot)  control  input 

Lateral  stick  driver  command  input 

Lateral  SAS  control  input 


Lateral  stick  trim  control  input 


aChf’aghr 


O'M 


oF 


Correction  to  "rotor  map"  derived  rotor 
drag  force  coefficient  arising  from 
angular  pitch  rate  at  rotor 

"Rotor  map"  power  coefficient  component 
due  to  compressibility  (reference 
value  for  Mach  number  at  which  maps 
were  calculated) 


"Rotor  map"  power  coefficient  component 
due  to  compressibility  (for  Mach 
number  at  which  aircraft  is  operating) 


* 


Mi 

In. 


In . 
In. 


In . 

In.  of 
equivalent 
cockpit 
control 

In  • 

In . 


In. 


In.  of 
equivalent 
cockpit 
control 

In. 
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Symbol 

A(D/q) 

Afe 

A (L/q) 
AP 

A(Y/q) 


Pesciiptian 

Nondimens ional  fuselage  drag  correction 
due  to  A^e 

Difference  between  actual  fuselage  fe 
and  fe  used  in  fuselage  aero  tables 

Nondimens ional  fuselage  lift  correction 
due  to  frfe 

Differential  SAS  port  pressure 

Nondimensional  fuselage  side  force 
correction  due  to  frfe 


Jlait 

Ft2 


Ft 

Ft2 

Lb/Ft^ 

Ft2 


A*QpiAfi^ 


B,0 


AQp  — 0 “ flp»  Afl^  ^ 

= Engine  speed  (-)  rotor  speed 

Second-order  upper  boost  actuator 
critical  damping  ratio 

Euler  pitch  angle  and  angular  rate 


Rad/Sec 


Rad 

Rad/Sec 


'.75 


eBF,0BR 


eCF’0CR 


Simulator  beam  pitch  angle 

Simulator  cockpit  pitch  angle 

.75  blade  radius  collective  pitch 
angle 

Longitudinal  cockpit  DCP  control  input 
expressed  in  terms  of  blade  pitch  angle 

Cockpit  collective  control  input 
expressed  in  terms  of  blade  pitch  angle 

Angle  between  load  X axis  and  helicopter 
X axis 


Deg 

Deg 

Rad 

Deg 

Deg 

Rad 


eRF,0RR 


eSF,0SR 


Cockpit  directional  control  input 
expressed  in  terms  of  blade  pitch  angle 

Cockpit  lateral  control  input  expressed 
in  terms  of  blade  pitch  angle 


Deg 

Deg 
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Symbol  Description 

8t„,©t_  Root  collective  pitch  at  full-down 

collective  stick 


Linear  rotor  blade  twist  (root  to  tip) 


eAP,0AS* 

0FP’eFS 


Combined  cockpit  and  SAS  mixed  controls 
at  aft  and  forward  swiveling  and 
pivoting  upper  boost  actuators 


eAPP,0ASP, 

0FPP*eFSP 


AP 


SAS, 


assas, 

0FP 


FS 


SAS, 

SAS, 


e 


of,  e 


oR 


Cockpit  controls  after  mixing  at  aft 
(or  forward)  rotor  pivoting  and 
swiveling  actuators 


SAS  controls  after  mixing  at  aft  (or 
forward)  rotor  pivoting  and  swiveling 
SAS  links 


Root  collective  pitch  (0^  is 
corrected  for  6 ^ ) 


0 r,w»9  nw  'Rotor  map"  .75  blade  radius 

oFM’  oRM  . . t . . . , 

collective  pitch 


oF' 


Actual  control  system  root  collective 
pitch  input 


e°F'6B 


Rate  of  change  of  rotor  collective 
pitch  with  longitudinal  stick  - forward 
rotor 


Rate  of  change  of  rotor  collective 
pitch  with  collective  stick  - forward 
rotor 


%5b 


Rate  of  change  of  rotor  collective 
pitch  with  longitudinal  stick  - aft  rotor 


Rate  of  change  of  rotor  collective 
pitch  with  collective  stick  - aft  rotor 


XF.XR 


Complete  inflow  ratio 


Unit 

Rad 

Rad 

Deg 

Deg 

Deg 

Rad 

Rad 

Rad 

Deg/In 

Deg/ In 

Deg/In 

Deg/In 
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Symbol 

Pesg.ription 

Unit 

>*V‘r 

Inflow  ratio  component  normal  to  the 
rotor  due  to  free-stream  velocity 

- 

xsl,xsl, 

XSL 

Lateral  cable  angle,  angular  velocity, 
and  anguJir  acceleration  (parallel  to 
YZ  plane  of  helicopter) 

Rad 

Rad/Sec 

Rad/Sec2 

Advance  ratio 

nR  hr 

- 

^FRM’^RRM 

"Rotor  map"  advance  ratio 

- 

^SL’^SL, 

Longitudinal  cable  angle,  angular 
velocity,  and  angular  acceleration 
(parallel  to  XZ  plane  of  helicopter) 

Rad 

Rad/Sec 

Rad/Sec2 

• 

y.SL’VSL’ 

V 

SL 

Lateral  differential  cable  angle, 
angular  velocity,  and  angular 
acceleration  (parallel  to  XY  plane 
of  helicopter) 

Rad 

Rad/Sec 

Rad/Sec2 

Po 

Air  density 

Slugs/Ft3 

P 

Air  density  at  sea  le\el 

Slugs/Ft3 

<J»CJp»CFR 

Rotor  solidity  ratio 

ttR 

- 

e(^aero 

Sum  of  forward  and  aft  rotor  aerodynamic 
torque  required 

Ft/Lb 

T 

Rotor  first  harmonic  flapping  time 
constant 

Sec 

Tx  ,TX 
XF  ar 

Rotor  inflow  lag  time  constant 

Sec 

Tn 

Engine  gas  generator  lag  time  constant 

Sec 

0,0 

Euler  roll  angle  and  angular  rate 

Rad 

Rad/Sec 

0c 

Simulator  cockpit  roll  angle 

Deg 
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Symbol 


Description 

Forward  rotor  6 3 control  phase  angle 

Special  control  phase  angle  used  to 
correct  "rotor  map"  cyclic,  pitch 
values  when  rotor  data  is  computed 
with  phase  angles  equal  to  the 

angle  instead  of  control  phase  angle 
Euler  yaw  angle  and  angular  rate 


+C 


n 


(qrf) 

(QRp) 


Simulator  beam  yaw  angle 

Simulator  cockpit  yaw  angle 

Direction  from  which  steady  wind 
blows  expressed  from  0°  to  £180' 

Undamped  natural  frequency  of  upper 
boost  actuator 

Engine  output  speed 

Reference  governor  speed 

Rotor  speed 

Tip  speed  of  the  forward  and  aft  rotors 


Unit 

Deg 

Deg 

Rad 

Rad/Sec 

Deg 

Deg 

Rad 

Rad/Sec 

Rad/Sec 

Rad/Sec 

Rad/Sec 

Ft/Sec 
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